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The  objectives  of  this  project  are: 

1.  Extend  the  receptivity  models  developed  for  2-D  disturbances  to  3-D  disturbances  to  treat 
receptivity  of  high-speed  boundary  layers  to  acoustics,  vorticity  and  entropy  fluctuations 
interacting  with  3-D  roughness/waviness. 

2.  Package  the  theoretical  model  in  a  computational  module  with  input  and  output  to  be 
accessible  for  a  standard  stability  solver  or  PSE  code. 

Analyses  of  high-speed  boundary  layer  receptivity  to  three-dimensional  external  waves 
interacting  with  roughness  element  or  waviness  were  perfonned.  The  FORTRAN  computational 
module  RECEPT.FOR  was  developed  and  delivered  to  the  end  user  Dr.  Roger  Kiminel  with  the 
usage  instructions.  This  module  allows  for  calculations  of  the  receptivity  coefficients  and  the 
weight  coefficients  of  unstable  modes.  It  was  shown  that  the  module  RECEPT  extracts  the 
nonnal  wave  of  relatively  small  amplitude  (-0.1%)  from  the  disturbance  field.  With  this  tool,  the 
user  can  identity  nonnal  waves,  which  latent  in  a  “messy”  disturbance  field  at  early  stages  of 
instability.  The  module  RECEPT  was  applied  to  the  problem  of  high-speed  boundary-layer 
receptivity  to  acoustic  disturbances  interacting  with  three-dimensional  roughness  element.  This 
option  was  demonstrated  by  numerical  examples.  Namely,  the  receptivity  coefficients  were 
calculated  for  two-dimensional  second  mode  and  three-dimensional  Tollmien-Schlichting  waves. 
Recommendations  for  design  of  receptivity  experiments  in  supersonic  quiet  wind  tunnels  were 
formulated.  This  helps  to  reduce  the  risk  and  cost  of  experiments,  especially  in  high-speed  wind 
tunnels  where  runs  are  quite  expensive. 
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Abstract 

Laminar-turbulent  transition  is  relevant  to  many  aspects  of  high-speed  vehicle  design.  Early 
transition  causes  a  significant  increase  in  heat  transfer.  Higher  heating  requires  a  higher- 
performance  thermal  protection  system  or  active  cooling  that  increases  cost  and  weight.  Vehicle 
drag,  engine  performance  and  vehicle  aerodynamics  are  also  affected.  Since  transition  prediction 
methods  currently  used  in  engineering  practice  are  subject  to  great  uncertainty,  conservative 
worst-case  conditions  are  normally  assumed  in  vehicle  design.  For  existing  or  near-term 
vehicles,  the  payoff  for  improved  transition  prediction  is  enhanced  perfonnance.  For  extended 
cruise  airbreathers,  transition  prediction  and  control  may  be  critical  to  achieving  design 
requirements.  Most  transition  prediction  methods  have  been,  at  best,  semi-empirical.  To  predict 
boundary-layer  transition  in  a  rational  way  and  capture  its  sensitivity  to  freestream  noise, 
turbulence,  surface  roughness,  steps/gaps/waviness,  etc.,  a  toolbox  combining  receptivity  theory 
with  linear  stability  solvers  is  needed.  This  project  addresses  this  need  with  theoretical  modeling 
of  high-speed  boundary-layer  receptivity  to  acoustic  disturbances  as  they  interact  with  3-D 
roughness  elements  and  waviness.  The  study  is  based  on  combined  asymptotic  and  numerical 
approach.  Using  the  multiple-modes  method,  the  receptivity  models  for  2-D  disturbances  are 
extended  to  the  3-D  case  to  handle  interaction  of  acoustics,  vorticity  and  entropy  spottiness  with 
3-D  roughness/waviness  of  aerodynamic  surfaces.  These  models  are  used  to  develop  a  compact 
and  computationally  robust  algorithm  providing  fast  calculations  of  receptivity  characteristics. 
The  numerical  algorithm  is  integrated  into  the  FOPRTRAN  module  RECEPT  that  can  be  linked 
with  a  standard  stability  solver.  This  module  allows  for  calculations  of  the  receptivity 
coefficients  and  weight  coefficients  of  various  modes  of  the  discrete  spectrum:  TS  waves,  second 
mode,  cross-flow  instability  etc. 
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Introduction 

Laminar-turbulent  transition  is  relevant  to  many  aspects  of  high-speed  vehicle  design  [1,2].  Early 
transition  causes  a  significant  increase  of  heat  transfer.  Higher  heating  requires  a  higher- 
performance  thermal  protection  system  (TPS)  or  active  cooling  that  increases  cost  and  weight. 
Even  modest  savings  in  TPS  weight  can  translate  into  significant  payload  increases  due  to  the 
low  payload  mass  fraction.  Vehicle  drag,  engine  performance  and  vehicle  aerodynamics  are  also 
affected.  Since  transition  prediction  methods  currently  used  in  engineering  practice  are  subject  to 
great  uncertainty,  conservative  worst-case  conditions  are  normally  assumed  in  vehicle  design. 
For  existing  or  near-term  vehicles,  the  payoff  for  improved  transition  prediction  is  enhanced 
performance.  For  extended  cruise  airbreathers,  transition  prediction  and  control  may  be  critical  to 
achieving  design  requirements. 

Boundary  layer  transition  is  difficult  to  predict  due  to  its  sensitivity  to  initial  conditions.  Most 
transition  prediction  methods  have  been,  at  best,  semi-empirical.  Recent  research  opens  a 
possibility  of  first-principles  transition  prediction  accounting  for  freestream  disturbance  levels 
and  surface  non-uniformities.  The  process  by  which  freestream  disturbances  interact  with  the 
boundary  layer  and  engender  instabilities  that  lead  to  turbulence  is  termed  receptivity.  This 
project  extends  past  basic  theoretical  research  on  receptivity  funded  by  EOARD,  AFOSR,  and 
AFRL  to  develop  the  theory  and  computational  modules  predicting  receptivity  of  the  boundary 
layer  on  a  rough  surface  exposed  to  acoustic  noise  or  other  external  disturbances  (vortical  waves 
and  entropy  spottiness). 

It  is  possible  to  identify  two  major  routes  to  turbulence.  The  first  route  is  observed  in  a  low- 
disturbance  environment  typical  for  flight  and  quiet  wind  tunnels.  In  this  case,  the  transition 
process  can  be  treated  as  an  initial  boundary-value  problem  that  comprises  excitation  of  unstable 
boundary-layer  modes  (receptivity  phase),  their  downstream  amplification  (linear  stability  phase) 
and  nonlinear  breakdown  to  turbulence  [3,  4].  The  second  route  is  observed  in  high-disturbance 
environments  and  associated  with  “bypass”  of  the  linear  phase  [5].  This  project  addresses  the 
first  route  of  transition  for  high-speed  boundary  layers. 

Reviews  on  transition  prediction  methods  [6-10]  indicate  good  progress  in  modeling  of  the  linear 
stability  phase  including  three-dimensional  stability  codes  under  the  quasi-parallel 
approximation  [11]  and  non-parallel  linear  stability  codes  based  on  the  multiple-scales  method 
[12-14],  Parabolized  stability  equations  (PSE)  methods  [15-17]  enable  capturing  of  both  non¬ 
parallel  and  weakly  non-linear  effects.  Direct  numerical  simulations  (DNS)  [18]  provide  insight 
into  breakdown  mechanisms.  Despite  this  progress,  the  operative  transition  prediction  method 
remains  the  eN  method  [10],  which  does  not  consider  receptivity  and  nonlinear  breakdown. 
These  phases  are  incorporated  into  the  value  N,  which  is  calibrated  for  a  particular  set  of 
experimental  data.  Malik  [10]  reported  that  the  eN  method  (with  N  pa  10)  provides  an  upper 
bound  on  the  transition  location  relevant  to  an  extremely  low-disturbance  limit  whereas  the 
lower  bound  (iV  «  1)  corresponds  to  bypass.  These  two  limits  are  often  wide  apart  and 
transition  occurs  anywhere  in  between.  To  predict  the  transition  locus  in  a  rational  way  and 
capture  its  sensitivity  to  freestream  noise  and  turbulence,  surface  roughness, 
steps/gaps/waviness,  etc.,  one  needs  a  more  sophisticated  toolbox,  which  combines  receptivity 
theory  with  stability  theory. 
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A  number  of  theoretical  models  have  been  developed  for  receptivity  of  high-speed  boundary 
layers  including  sharp  leading-edge  receptivity  to  freestream  acoustic  and  vortical  disturbances 
[19,  20],  receptivity  to  vortical/entropy  disturbances  near  the  lower  branch  of  the  neutral  curve 
[21],  local  and  distributed  receptivity  to  wall-induced  disturbances  such  as  skin  vibrations 

[22.23] ,  and  receptivity  to  external  disturbances  interacting  with  a  small  roughness/waviness 

[22.24] .  The  multiple -modes  method  described  in  [24]  (see  also  [21])  serves  as  a  general 
formalism  for  the  aforementioned  theoretical  models.  This  provides  a  good  launching  pad  for 
further  extension  of  receptivity  studies  to  more  practical  cases  accessible  for  experimental 
verification. 

Important  mechanism  of  boundary-layer  receptivity  is  the  interaction  of  acoustic  disturbances 
with  steady  disturbances  induced  by  spatial  variations  of  boundary  conditions  on  the  wall  surface 
such  as  roughness/waviness,  local  heating/cooling  and  suction/blowing.  Early  experimental  [25] 
and  theoretical  [26-28]  studies  of  this  mechanism  were  performed  in  Russia  for  subsonic  flows. 
Their  results  were  summarized  in  Ref.  [24].  An  asymptotic  theory  of  the  Tollmien-Schlichting 
(TS)  wave  excitation  due  to  interaction  of  acoustic  waves  with  surface  irregularities  was 
developed  in  Refs.  [29-31]  and  reviewed  in  Refs.  [31-33].  Detailed  parametrical  studies  were 
discussed  in  Refs.  [34-37],  The  theoretical  predictions  are  in  good  agreement  with  the 
experiment  [38]  conducted  on  a  flat  plate  in  low-speed  wind  tunnel. 

Choudhari  and  Streett  [39]  analyzed  receptivity  of  high-speed  boundary  layers  (in  the  Mach 
number  range  from  1.15  to  4.5)  to  freestream  disturbances  interacting  with  surface  irregularities. 
They  reported  that,  although  the  Mach=4.5  boundary  layer  is  nearly  6  times  thicker  than  in  the 
Mach=0  case,  the  instability  amplitudes  generated  by  a  wall  hump  of  fixed  height  are  only 
marginally  smaller  than  those  at  low  speeds.  This  example  indicates  that  receptivity  of  high¬ 
speed  boundary  layer  to  acoustic  plus  roughness  disturbances  may  play  important  role  in  the 
initial  phase  of  transition. 

In  many  cases,  it  is  needed  to  identify  a  particular  mode  of  discrete  spectrum  in  a  disturbance 
field,  which  is  resulted  from  DNS  or  experimental  measurements.  This  is  not  a  trivial  problem 
because  the  disturbance  field  may  contain  many  modes  of  discrete  and  continuous  spectra.  The 
unstable  mode  may  be  latent  in  such  a  “messy”  field,  especially  in  early  phases  of  instability 
evolution  when  the  amplitude  of  unstable  waves  is  relatively  small.  This  motivated  us  to  include 
in  the  receptivity  module  an  option  providing  calculations  of  the  weight  coefficients  for 
disturbances  of  discrete  spectrum  (TS  waves,  second  mode  etc.) 

In  this  report,  we  describe  the  FORTRAN  receptivity  module  RECEPT.FOR  and  discuss  basic 
features  of  receptivity  to  acoustic  waves  interacting  with  roughness  elements.  In  Section  1,  we 
outline  basic  relations  incorporated  into  the  module.  The  theoretical  analysis  leading  to  these 
relations  is  presented  in  Appendix.  In  Section  2,  we  discuss  the  module  input  and  output  and 
explain  how  to  couple  the  module  with  a  linear  stability  solver.  In  Section  3,  we  illustrate  how  to 
use  the  module  output  for  calculations  of  physical  quantities.  In  Section  4,  we  present  numerical 
examples  demonstrating  computations  of  the  weight  coefficients.  In  Section  5,  we  consider 
receptivity  for  2-D  second  mode  generated  by  3-D  acoustic  waves  interacting  with  local  3-D 
roughness  element.  For  Mach=6  boundary  layer  on  a  flat  plate,  we  present  numerical  results 
illustrating  distributions  of  the  receptivity  coefficient  versus  the  angles  of  inclination  and 
incidence  of  fast  and  slow  acoustic  waves.  In  Section  6,  we  discuss  receptivity  for  3-D  Tollmien- 
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Schlichting  (TS)  waves.  Using  the  module  RECEPT  we  calculate  the  receptivity  coefficient  for 
various  angles  of  acoustic  and  TS  waves.  Special  attention  is  paid  to  the  fact  that  the  receptivity 
coefficient  is  singular  if  the  resonant  roughness  component  is  transonic;  i.e.,  the  freestream 
velocity  normal  to  the  waviness  front  equals  to  the  speed  of  sound.  In  Section  7,  we  fonnulate 
recommendations  for  design  of  receptivity  experiments  in  quiet  high-speed  wind  tunnels.  In 
Section  8,  we  conclude  the  report,  and  in  Section  9  we  discuss  the  future  effort. 


1.  Basic  relations 

Herein  we  briefly  summarize  basic  analytical  relations,  which  are  incorporated  into  the 
FORTRAN  module  RECEPT. FOR.  Details  of  the  theoretical  analysis  are  given  in  Appendix. 

Consider  the  boundary  layer  on  a  2-D  body  (planar  or  axisymmetric),  on  a  swept  wing  or  similar 
configurations  for  which  the  boundary  layer  is  predominantly  two-dimensional.  At  a  certain 
streamwise  station  x  =  x* ,  the  mean-flow  profiles  are 

U  =  U*/U;=U(y),W  =  W*/U:=W(y),  (1.1) 

T  =  r/T;  =  T(y). 


Hereafter  the  longitudinal  x* ,  normal  y* ,  transversal  z*  coordinates  and  time  t*  are  made 
nondimensional  using  the  boundary-layer  length  scale  6*  and  the  velocity  U*  as 


(x,y,z)  =  (x*,  y*,  z*)/8*,t 


t*U*  /  8* 


(1.2) 


where  subscript  “e”  denotes  quantities  at  the  upper  boundary-layer  edge.  Introduce  the  small 
parameter  e  =  b*  /  x* ,  which  is  a  ratio  of  the  boundary-layer  scale  to  the  global  scale  x*0 .  If  the 

boundary-layer  scale  is  specified  as  6*  =  /  U*  ,  then  the  Reynolds  number  is 

R  =  \]u*ex*Q  /  v*  =  e~l .  The  slow  streamwise  variable  is  defined  as  x1  =  ex .  The  vertical 
velocity  component  V  =  V*  /  U*  =  0(e)  is  neglected  in  the  first-order  approximation. 

Disturbances  of  the  velocity  components  (■ u,v,w ),  pressure  p  and  temperature  9  are 
nondimensionalized  as 


(u,v,w) 


/  *  *  *  \  * 

(u  u  ,w  )  „  _  1  n 

ue  P;uf 


(1.3) 


1.1  Disturbance  field 

A  disturbance  of  angular  frequency  u  is  given  by  the  vector  function 
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oy  oy  oy 


(1.4) 


The  disturbance  field  comprises:  the  stationary  disturbance  induced  by  the  wall  roughness  of  the 
amplitude  h ;  the  external  disturbance  of  continuous  spectrum  with  the  amplitude  sa ;  modes  of 
discrete  spectrum 

*P  =  hQw(x,y,z )  +  £0Fa (xvy)  exp (iaax  +  if3az  -  iut )  +  F(x,y,z)exp(-iut)  +  (c.c) .  (1.5) 


The  disturbance  of  continuous  spectrum  may  be  acoustic  wave,  vortical  wave  or  entropy  wave. 
Hereafter  these  waves  are  called  “external  waves”.  Modes  of  the  discrete  spectrum  are  expressed 
in  the  form 


F  =  h£aJ2lFno(XVy)  +  £Fnl(XVy )  +  -]eXP(*£  X  +  */?*)  >  (  1 -6) 

n 

x i 

Sn  =  f  o:n(x1,cj,/3)dx1 ,  an  =  an(xvu,P),  (1.7) 


where  the  amplitude  of  Mode  is 

F0  =  c^xjA^MaJ,  an  =  an(xvu,  (3) . 


(1.8) 


1.2  Resonant  excitation  of  a  discrete  mode  by  an  external  wave  interacting 
with  the  wall  waviness 


The  stationary  disturbance  induced  by  a  wavy  wall  is  expressed  as 

Q»  =  Fw(xvy)exp(iawx  +  if3wz)  +  (c.c) . 


(1.9) 


At  the  station  x  =  x  ,  the  roughness  wave-number  is  determined  from  the  resonance  conditions 


=  aB(s10, <*>,/?)  -  aa ,  (3W  =  (3  -  f3a . 


The  amplitude  of  Mode  “n”,  propagating  downstream  from  the  resonance  point  x  ,  is 


(1.10) 


F „  =  £ahcn  (xl0,  P)An(xv  y)  exp 


1 

zs-1  J  an(x1)dx1  +  ifiz 


xi  >  ®10  ’ 


Cn(XW’P)  =  £  ? 


2m 


(1.11) 

(1.12) 
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,  9an 
b  =  —JL(xJ. 


dx, 


q  = 


i<B0,G(xJF„> 

OH 


10’  o/o 

«0  =  an(xw,uJ,P),  A  =  A„(x w,y,a  ),  B()  =  B„(a:  10,y,a) . 


(1.13) 

(1.14) 

(1.15) 


Here  cn  is  receptivity  coefficient,  which  couples  the  amplitudes  of  external  disturbance  and 
wavy  wall  with  the  amplitude  of  Mode  n  at  the  station  x1  =  xw .  Derivation  of  these  relations  and 
definitions  of  other  quantities  are  given  in  Appendix. 

1.3  Excitation  of  a  discrete  mode  by  an  external  wave  interacting  with  a 
local  hump 

It  is  assumed  that  the  surface  roughness  is  a  local  hump  of  the  shape 

yw(x,y,z)  =  hg{x',z'),  x  =  X-X0,  z  =  z-zQ.  (1.16) 

Fourier  transform  of  g(x',  z')  is  expressed  as 


+00  +0O 

p(aw,Pw)  =  f  dz'J  dx'g{x',z')exp{-iawx'  -  i(3wz') . 

— OO  — OO 


(1.17) 


For  a  two-dimensional  hump  with  g  =  g{x')  and  /3W  =  0  ,  Mode  n  has  the  amplitude 


K  =  ^ahcn(xm,p)An(x/,y)exp 


i  J  an{x')dx'  +  i/3z 


xi>xw: 


cn(xw,(3)  =  qp(aj, 

aw0  =  an  Oi0>  P)  ~  aa  »  P  =  Pa  • 


For  a  three-dimensional  hump  with  g  =  g(x',z').  Mode  n  has  the  amplitude 


(1.18) 

(1.19) 

(1.20) 


_  7  -t-OO 

F"  =  2nf  Cn(X10’fi)An(Xl’y;fi) 


exp 


1 

ie _1  J  an(x1)dx1  +i/3z' 


d/3,  x1>  xw ; 


cn(x10>P)  =  VPip^PJ’ 
o  =  an(X  10>V>P)  -  aa  ,  P  =  Pw  +  Pa 


(1.21) 

(1.22) 

(1.23) 
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1.4  Weight  coefficients 

In  many  cases,  it  is  needed  to  identify  a  mode  of  the  discrete  spectrum,  say  Mode  n,  in  a 
disturbance  field,  which  is  resulted  from  DNS  or  experimental  measurements.  This  is  not  a  trivial 
problem  because  the  disturbance  field  may  contain  many  modes  of  discrete  and  continuous 
spectra.  The  amplitude  of  unstable  mode  may  be  latent  in  such  a  “messy”  field,  especially  in 
early  stages  of  instability. 

The  weight  coefficient  of  a  particular  discrete  mode  is  calculated  as  follows.  The  disturbance 
field  is  assumed  to  be  known.  Performing  Fourier  transform  with  respect  to  t  and  z ,  we  obtain 

Q  =  F(x,  y)exp(i/3z  —  iut) .  (1-24) 

This  Fourier  component  may  contain  modes  of  discrete  spectrum  (TS  wave,  second  mode,  etc.), 
which  are  solutions  of  the  eigenvalue  problem.  If  the  disturbance  field  at  a  particular  station 
x  =  xQ  is  F  (y)  =  F (x  ,y) ,  then  Mode  n  of  the  eigenvalue  an(xQ,u ;,/?)  has  the  amplitude 


F„  =  cn(x0,p)An(x0,y,an), 

where  the  weight  coefficient  is  expressed  as 


cn(®0,/?) 


< 


0 

da 


< 


d_H o 

da 


(®o > 


(1.25) 


(1.26) 


1.5  Summary 

Basic  relations  of  Sections  1.1  and  1.2  have  similar  structures.  This  allows  for  the  development 
of  a  compact  computational  module,  which  calculates: 

•  Receptivity  coefficients  cn{x  ,0)  given  by  (1.12),  (1.22)  for  the  case  of  receptivity  to 

external  wave  interacting  with  roughness-induced  disturbance 

•  Weight  coefficient  cn(x  ,0)  given  by  ( 1 .26). 

This  module  can  be  linked  with  a  standard  linear  stability  solver  as  described  in  Section  2. 


2.  Receptivity  module  RECEPT  and  usage  instructions 

The  FORTRAN  module  RECEPT. FOR  calculates: 

•  The  coefficient  (1 . 12)  for  receptivity  to  external  wave  interacting  with  a  wavy  surface 

•  The  coefficient  (1.22)  for  receptivity  to  external  wave  interacting  with  a  local  hump 

•  The  weight  coefficient  ( 1 .26)  for  a  given  disturbance  field 
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In  this  version,  it  is  assumed  that  the  branch  points  of  the  discrete  spectrum  are  far  from  the 

OH 

resonance  region;  i.e.  <  0  (see  Appendix). 

da 


The  RECEPT.FOR  is  a  FORTRAN  subroutine,  which  needs  to  be  linked  with  a  linear  stability 
solver.  In  order  to  specify  the  full  input  for  this  module,  the  user  should  follow  steps  indicated  in 
Sections  2. 1-2.3.  Note  that  not  all  of  the  input  data  are  needed  for  a  particular  option.  Subsets  of 
data  for  different  options  are  given  in  Section  2.4. 

2.1  Mean  flow 

2.1.1  Specify  gas  parameters: 


PR  -  Prandtl  number  (real) 

G=  7  -  specific  heat  ratio  (real) 
3  (bulk  viscosity) 


K=  e  = 


(real) 


2  (viscosity) 

(K=0.  corresponds  to  Stokes  hypothesis;  K=1.2  was  recommended  by  L.  Mack) 


2.1.2  Specify  the  mean-flow  profiles  and  their  derivatives  with  respect  to  y  at  station  xl  =  x 


du  dw  dT 

(■ U,W,T ),  {DU.DW ,DT)  = 

ay  dy  dy 
dy  dy  dy 


(2.1) 


as  well  as  the  distributions  of  viscosity  and  its  derivatives  with  respect  to  the  gas  temperature 


a_ 

* 

M e 


M  >  v>\y) 


~dT 


(: V ) >  Av) 


(2.2) 


2.1.3  Specify  the  local  Mach  number  M  =  — %-  at  the  upper  boundary  y  =  yr  (real) 

ae 

2.1.4  Store  these  data  in  the  common  blocks: 

COMMON/OUT/K,PR,G,M 

C0MM0N/BASA5/NBAS,YY,BU,BU1,BU2,BT,BT1,BT2 
*,BW,BW  1  ,BW2,BMU,BMU  1  ,BMU2 


where 

YY(I)=y  (real),  1=1,..., NBAS; 
BU(I)=t/  (real),  1=1,..., NBAS; 
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BU1(I )=DU  (real),  1=1,..., NBAS; 

BU2(I)=  D2U  (real),  I=1,...,NBAS; 

BT(I)=T  (real),  1=1,..., NBAS; 

BT1(I )=DT  (real),  1=1,..., NBAS; 

BT2(I)=  D2T  (real),  1=1,..., NBAS; 

BW(I)=1R  (real),  1=1,..., NBAS; 

B W 1  (I)=  D  W  (real),  1=  1 , . . .  ,NB AS ; 

BW2(I)=  D2W  (real),  1=1,..., NBAS; 

BMU(I)=  ii  (real),  1=1,...  ,NB AS ; 

BMU 1  (I)=  (real),  1=1,...  ,NB AS ; 

dT 

t2 

BMU2(I)=  — ^  (real),  1=1,..., NBAS; 
d  T3 

NBAS  =number  of  grid  points  in  the  range  0  <  y  <  ye ;  ye  is  upper  boundary  of  the 
integration  domain. 

2.2  Solutions  of  direct  and  adjoint  problems 

2.2.1  Specify  the  disturbance  angular  frequency  u> ,  the  wave-number  component  (3  and  the 
Reynolds  number  R .  Solve  the  direct  eigenvalue  problem  for  a  particular  mode  of  discrete 
spectrum  (say  Mode  “«”)  at  the  station  x1  =  x 

dA 

=  H0 A  ,  (2.3) 

dy 

A  =  A  =  A  =  at  =  0’  y  =  °> 

A  — >  0,  y  — >  oo  . 

Store  the  eigenvalue  a  =A,  and  the  parameters  u  =W,  [3  =B  and  R  in  the  common  block 
COMMON/SCAL2/A,B,W,R 
where  A,B,W,R  are  complex. 

2.2.2  Calculate  the  eigenvector  A (y)  of  the  direct  problem  (2.3)  in  the  range  0  <y<ye. 

2.2.3  Calculate  the  asymptotic  form  of  A (y)  for  y  >  ye 

A  =  Ao?)  exP [Aj(2/  -  Ve)\>  (2-4) 

3= 1 

where  Re(A  )  <  0 . 
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2.2.4  Solve  the  adjoint  problem  using  a ,  (3  and  u>  of  the  direct  problem  (2.3) 


dB 

dy 


o,  y 


B 


0,  y  =  0, 


(2.5) 


2.2.5  Calculate  the  eigenvector  B  in  the  range  0  <y<ye. 

2.2.6  Calculate  the  asymptotic  form  of  B  for  y  >  ye 

B  =  J2Bo]  exP -  Ve)\  •  (2-6) 

I 

2.2.7  Store  solutions  of  the  direct  and  adjoint  problems  in  the  common  block 

COMMON/SCAL 1/  YZ,Z,DZ,ZO,DZO,LZ,LDZ,SC,NZ,MF,L 

where 


YZ(K )=y  (real),  K=1,...,NZ,  is  vector  of  grid  points; 

Z(I,K)=A  (complex),  1=1,.  ..,8  is  component  number,  K=1,...,NZ; 

DZ(I,K)=B  (complex),  1=1,. ..,8  is  component  number,  K=1,...,NZ; 

Z0(I,J)=  Aj^  (complex),  1=1,.  ..8,  J=  1 , ... ,4  ; 

DZ0(I,J)=  Bjj’  (complex),  1=1,... 8,  J=l,... ,4  ; 

LZ=LDZ=A/  (complex),  J=l,..., 4; 

NZ  -  (integer)  is  number  of  grid  points  in  the  range  0  <  y  <  ye  for  the  eigenvectors  Z  and  DZ; 


SC  (complex),  MF  (integer)  and  L  (integer)  are  additional  parameters,  which  are  not  specified  by 
the  user. 


2.2.8  Calculate  the  derivative  DADX1= 


da 

dx1 


This  can  be  done  by  solving  the  direct 


problem  (2.3)  at  two  close  stations  x1  =  x  —  Ax1  and  x1  =  x  +  Axl .  Then 


da 

dx± 


a(xw  +  Ax1)-a{xw  -  AxJ 
2Axx 


(2.7) 


The  parameter  DADX1  (complex)  is  a  formal  parameter  of  RECEPT. 
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2.3  Solution  for  external  wave 

Hereafter,  we  discuss  the  case  when  the  external  disturbance  is  an  acoustic  wave.  The 
replacement  of  acoustic  wave  by  vortical  or  entropy  wave  is  straightforward. 

2.3.1  Solve  the  local  problem  for  the  acoustic  wave  of  the  wave  numbers  aa  ,  (3a  and  frequency 
u  (this  frequency=frequency  of  the  direct  and  adjoint  problems  of  Section  2.2)  at  the  station 


^Aa  —HA 

dy 

4,  ='4.«  =  4S  =  AI  =  °.  »  =  °. 

external  conditions  at  the  upper  boundary  y  =  ye. 


(2.8) 


The  external  condition  may  be  different.  For  example:  (the  pressure  amplitude  of  incident 
acoustic  wave  at  y  =  ye  )=1 . 

2.3.2  Calculate  the  external-disturbance  vector  A a(y)  in  the  range  0  <  y  <  ye 

2.3.3  Calculate  the  asymptotic  form  of  A a{y)  for  y  >  ye 

Aa  =  eXP (KjiV  -Ve))>  (2-9) 

3= 1 

where  R,e(A.)  <  0.  Note  that  the  external-disturbance  vector  is  a  sum  of  five  vectors.  Three  of 

them  exponentially  decay  as  y —>  oo.  Two  others,  with  Re(Aa)  =  0,  oscillate  outside  the 
boundary  layer  without  damping.  They  represent  an  external  disturbance  outside  the  boundary 
layer.  For  the  acoustic  disturbance,  the  oscillatory  components  represent  the  incident  and 
reflected  waves. 


2.3.4  Store  characteristics  of  the  external  disturbance  in  the  common  block 


COMMON/SCALA/  YA,ZA,ZAO,LA,AAK,BAK,NA 


where 

YA(K )=y  (real),  K=1,...,NA,  is  vector  of  grid  points; 

ZA(I,K)=  Aft  (complex),  1=1,.  ..,8  is  component  number,  K=1,...,NA; 
ZA0(I,J)=A9)  (complex),  1=1,... 8  is  component  number,  J=l,...,5; 

LA=A/  (complex),  J=l,..., 5; 

AAK=oia  (complex); 

BAK=/30  (complex); 

NA  -  (integer)  is  number  of  grid  points  in  the  range  0  <  y  <  ye  for  the  vector  ZA. 
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2.4.  Solution  for  stationary  disturbance  induced  by  wavy  wall 

2.4.1  Solve  the  problem  for  the  stationary  disturbance  induced  by  the  wavy  wall  of 

aw  =  a  -  aa,  (3W  =  (3  -  (3a  and  uw  =  0  at  the  station  x1  =  xw 


~Tl  =  H  0Vw,  (2.10) 

dy 

F  =  —D  U,  F  =  - DT ,  F  _=-DW,y  =  0  , 

F»  -»•  0,  y  -»•  oo . 

2.4.2  Calculate  the  vector  function  F w(y)  in  the  range  0  <y<ye. 

2.4.3  Calculate  the  asymptotic  form  of  F w(y)  for  y  >  ye 

»-».)]■  <2-n> 

i=i 

where  Re(A  )  <  0  . 

v  WJ J 

2.4.4  Store  characteristics  of  the  stationary  disturbance  in  the  common  block 

COMMON/SCALV/  YV,ZV,ZVO,LV,AVV,BVV,NV 


where 

YV(K)=  y  (real),  K=1 , . . .  ,NV,  is  vector  of  grid  points; 

ZV (I,K)=  F(i  (complex),  1=  1 , ... ,8  is  component  number,  K=  1 , . . .  ,NV ; 

ZVO(I,J)=F'(0)  (complex),  1=1,. ..8  is  component  number,  J=l,...,4; 

LV=  A  (complex),  J=  1 , . . . 4; 

A V V=  aw  (complex); 

BVV=  j3w  (complex); 

NV  -  (integer)  is  number  of  grid  points  in  the  range  0  <  y  <  ye  for  the  vector  ZV. 

2.5  Disturbance  field 

2.5.1  Specify  the  disturbance  field  vector  Ffj(i/)  in  the  range  0  <  y  <  ye  at  the  station  x1  =  x1Q. 

2.5.2  Store  the  vector  F0(y)  in  the  common  block 

COMMON/FIELD/YAO,AO,NAO 


where 
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YAO(K )=y  (real),  K=1,...,NA0,  is  vector  of  grid  points; 

A0(I,K)=  FQ  (complex),  1=1 , ...  ,8  is  component  number;  K=1 , . . .  ,NA0; 

NAO  -  (integer)  is  number  of  grid  points  for  the  vector  AO. 

2.6  Options  of  RECEPT.FOR 

The  module  includes  the  subroutine  RECEPT  and  the  auxiliary  subroutines  and  functions: 

•  Subroutine  SCAL  (calculations  of  scalar  products) 

•  Subroutine  BASIC  (interpolation  of  mean-flow  profiles) 

•  Function  PARAB2  (interpolation  of  real  function) 

•  Function  PARAB3  (interpolation  of  complex  function) 

•  Function  SIM  (integration) 

Calculations  are  performed  with  double  precision: 

•  All  real  quantities  are  described  as  REAL  *8 

•  All  complex  quantities  are  describes  as  COMPLEX*  16 

The  calling  line  is 

CALL  RECEPT (ICONTRJGNOR, Y 1  ,N,DADX  1  ,EPSIL,CC) 


Since  the  major  portion  of  input  data  are  specified  in  COMMON  blocks,  the  parameter  NNY  is 
included  in  RECEPT  and  BASIC  to  reserve  space  for  vectors.  The  user  needs  to  specify  NNY  in 
these  subroutines:  NNY  should  be  larger  then  NA,  NZ,  NV  and  NBAS. 


2.6.1  To  calculate  the  receptivity  coefficient  (1.12)  for  external  waves  interacting  with  the  wavy 
wall,  the  user  should  put  ICONTR=l  and  specify  the  input  data  of: 


Steps  1.1 -1.4  (mean-flow  data) 

Steps  2. 1-2.8  (solutions  of  the  direct  and  adjoint  problems  for  particular  mode  of  discrete 

dcY 

spectrum,  as  well  as  the  derivative  DADX1= - (x  )  (complex)) 

dx. 


•  EPSIL=  e  (real) 

•  Steps  3 . 1  -3 .4  (external  wave) 

•  Steps  4. 1-4.4  (stationary  disturbance  from  wavy  wall) 

•  Yl= ye  -  upper  boundary  (real) 

•  N  -  number  of  knots  for  calculation  of  the  scalar  product  integrals  in  the  range  0  <  y  <  ye 
(integer,  odd) 


22 


ISTC  2393 


If  the  linear  stability  solver  does  not  provide  the  asymptotic  forms  (2.4),  (2.6)  and  (2.9),  (2.11), 
then  the  user  should  put  IGNOR=l.  In  this  case,  the  contribution  to  the  scalar  products  from  the 
region  ye  <  y  <  oo  is  ignored.  If  the  upper  boundary  y  =  ye  is  not  large  enough,  this  may  lead 
to  significant  error  in  calculations  of  the  receptivity  coefficients.  To  evaluate  this  error  the  user 
should  perform  calculations  at  different  ye . 

If  the  asymptotic  fonns  (2.4),  (2.6)  and  (2.9),  (2.11)  are  available,  then  the  user  should  put 
IGNOR=0.  In  this  case,  it  is  recommended  to  specify: 

•  The  upper  boundary  as  Yl=r/e  «  3A,  where  A  is  nondimensional  displacement 
thickness 

•  The  number  of  knots  N=5  01 

However,  this  recommendation  is  not  universal.  For  different  boundary-layer  modes  and 
external  waves,  the  optimal  ye  may  be  different.  The  user  needs  to  check  sensitivity  of  the 
computational  result  to  Y 1  and  N. 

Output  of  this  option  is 

CC  -  receptivity  coefficient  (1.12)  (complex) 

2.6.2  To  calculate  the  receptivity  coefficient  (1.22)  for  external  waves  interacting  with  a  local 
hump,  the  user  should  put  ICONTR=2  and  specify  the  input  data  of: 

•  Steps  1.1 -1.4  (mean-flow  data) 

•  Steps  2. 1-2.7  (solutions  of  the  direct  and  adjoint  problems  for  particular  mode  of  discrete 
spectrum) 

•  Steps  3 . 1  -3 .4  (external  wave) 

•  Steps  4. 1-4.4  (stationary  disturbance  from  wavy  wall) 

•  Y1 =ye  -  upper  boundary  (real) 

•  N  -  number  of  knots  for  calculation  of  the  scalar  product  integrals  in  the  range  0  <  y  <  ye 
(integer,  odd) 

In  this  case,  the  formal  parameter  DADX1  and  EPSIL  are  not  used. 

If  the  LST  solver  does  not  provide  the  asymptotic  forms  (2.4),  (2.6),  (2.9),  (2.11),  then  the  user 
should  put  IGNOR=l  and  follow  instructions  of  Item  2.6.1. 

Output  of  this  option  is 

CC  -  receptivity  coefficient  ( 1 .22)  for  external  wave  plus  local  hump  (complex) 
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2.6.3  To  calculate  the  weight  coefficient  of  the  discrete  mode  in  the  disturbance  field,  the  user 
should  put  ICONTR=3  and  specify  the  input  data  of: 

•  Steps  1.1 -1.4  (mean-flow  data) 

•  Steps  2. 1-2.7  (solutions  of  the  direct  and  adjoint  LST  problems) 

•  Steps  5 . 1  -5 .2  (disturbance  field) 

•  Y1 =ye  -  upper  boundary  (real) 

•  N  -  number  of  knots  for  calculation  of  the  scalar  product  integrals  in  the  range  0  <  y  <  ye 
(integer,  odd) 

In  this  option,  the  asymptotic  fonns  (2.6),  (2.8)  and  the  parameters  DADX1,  IGNOR  are  not 
used.  Contributions  to  the  scalar  products  from  the  region  ye  <  y  <  oo  are  neglected.  If  the 
upper  boundary  y  =  ye  is  not  large  enough,  this  may  lead  to  significant  error  in  calculations  of 
the  weight  coefficient.  With  this  option,  the  user  should  perform  calculations  at  different  ye  to 
evaluate  the  error  (see  Section  4). 

Output  of  this  option  is 

CC  -  the  weight  coefficient  ( 1 .26)  (complex) 


3.  How  to  use  the  module  output? 

Using  the  receptivity  coefficient  CC  and  the  eigenfunction  vector  Z(I,J)  of  discrete  mode,  the 
user  can  calculate  the  amplitude  of  any  physical  quantity.  For  example,  the  pressure  disturbance 
on  the  wall  p(y  =  0)  at  station  x  =  xQ  is  calculated  as: 


For  wavy  wall:  p  =  hea  Re[CC  x  Z( 4,1)]  +  (c.c) 


For  local  hump:  p  =  hea  Re 


-J  +oo 

—  f  CC(/3)  x  Z(4,l)  x  p(a  ,/3)  x  exp(i/3z)d/3 
2ir  J 


(3.1) 

(3.2) 


where  1=4  corresponds  to  the  pressure  component  of  the  vector  Z(I,J)  (see  (1.4)),  and  J=1 
corresponds  to  y  =  0  . 

For  analysis  of  hot-wire  measurements,  the  user  needs  to  calculate  the  mass-flow  disturbance, 
which  is  expressed  as 


Qm(y) 


u{y)  + 


7  M'2p{y) 


%) 


Ay) 


U(y ) 


Ay) 


(3.3) 


At  the  point  y  =  y  the  correspondent  quantity  is  calculated  using  the  components  of  vector  Z  as 
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QM(J)=(Z(1,J)+(G*M*M*Z(4,J)-Z(5,J)/T(J))*U(J))/T(J).  (3.4) 

Replacing  Z(4,l)  by  QM(J)  in  (3.1)  and  (3.2),  the  user  obtain  the  mass-flow  amplitude  of  the 
discrete  mode. 

Similar  relations  are  used  for  calculations  of  the  discrete  mode  amplitude  in  the  disturbance  field 
(the  option  ICONTR=3). 


4.  Weight  coefficients 

The  discrete-mode  weight  coefficient  is  calculated  using  RECEPT  with  the  option  ICONTR=3. 
Herein  this  option  is  demonstrated  by  numerical  examples. 

Consider  the  mean  flow  in  the  boundary  layer  on  a  flat  plate.  This  flow  is  modeled  using  the  self¬ 
similar  (compressible  Blasius)  solution,  which  provides  the  longitudinal  velocity  profile 
U(y)  =  U*  /U*  and  the  temperature  profile  T(y)  =  T*  /  T* ,  where  the  vertical  coordinate  is 

v  =  v'  /  PP  1  K  ■  The  mean-flow  pressure  is  constant  across  the  boundary  layer, 
P  =  P*  /  P*  =  1 .  The  fluid  is  a  perfect  gas  with  Prandtl  number  Pr=0.72  and  the  specific  heat 
ratio  7  =  1.4  ;  viscosity  depends  on  temperature  as  //  =  //*  (T*  /  T'f" ;  the  second  viscosity  is 
zero.  The  local  Mach  number  is  Me  =  6  ;  the  plate-surface  temperature  equals  to  the  temperature 
of  adiabatic  wall,  T*  =  T*  .  The  profiles  U(y)  and  T(y)  are  shown  in  Figs.  la,b.  The  boundary- 
layer  edge  (at  which  U  =  0.99  )  is  approximately  y  =  16  (see  the  red  line  in  Fig.  la). 

The  disturbance  field  F '  (y)  of  the  frequency  parameter  F  =  1.3  x  10  4  is  specified  at  the  fixed 

station  x* ,  which  corresponds  to  the  Reynolds  number  R  =  ^U*x*  /  v*e  =  103 .  This  field 
contains  the  acoustic  wave  of  amplitude  F a(y)  and  the  second  mode  of  amplitude  c2F2(y),  where 
c2  is  the  weight  coefficient;  i.e.,  the  disturbance  field  distribution  is  specified  as 

F,W  =  F,(!/)  +  c!F2(!/).  (4.1) 

The  acoustic  component  F a(y)  represents  2-D  fast  acoustic  wave  of  the  angle  of  incidence 

0  =  60° .  Its  phase  velocity  is  Ca  =  —  =  1  4 - - - .  The  second-mode  component  has 

aa  M.  cost-) 

the  eigenvalue  a2  =  1.38896E-01-/9.26306E-04,  which  corresponds  to  2-D  unstable  nonnal 
wave  in  the  boundary  layer. 

In  our  calculations,  we: 

•  Assign  c2  and  calculate  the  disturbance  field  (4.1) 

•  Call  RECEPT  with  the  option  ICONTR=3  and  calculate  the  weight  coefficient  c2  mm 
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•  Compare  the  weight  coefficient  c2num  predicted  by  RECEPT  with  the  actual  weight 
coefficient  c2 

The  calculations  are  conducted  using  a  uniform  grid  versus  y  with  541  grid  points.  The  upper 
boundary  yi  =  ye  of  the  computational  domain  0  <y<y1  is  varied  in  order  to  investigate  the 
sensitivity  of  c„  to  y, . 

Figure  2  shows  the  pressure  distributions  pr  (y)  of  the  second  mode  with  c2  =  1  (red  line),  for 

the  acoustic  wave  (blue  line),  and  for  the  disturbance  field=(acoustic  wave)+(second  mode) 
(black  line).  This  example  corresponds  to  the  upper  boundary  y  «  42 .  The  weigh  coefficient 

predicted  by  the  module  RECEPT  is  c  =  1.00916.  With  this  coefficient,  the  second-mode 

eigenfunction  extracted  from  the  disturbance  field  (red  circles)  practically  coincides  with  the 
actual  second-mode  component  (red  line). 

Figure  3  shows  the  numerical  weigh  coefficient  c2  as  a  function  of  the  actual  weight 
coefficient  c„  for  various  upper  boundaries  y, .  As  c„  decreases,  the  curve  u  (c)  deviates 
from  the  bisectrix  c  =  c2  that  indicates  the  computational  error.  For  the  upper  boundary 
y  =  30  2  A  (where  A =16  is  the  boundary-layer  thickness),  the  error  is  ~  2.  x  10  2 .  It  quickly 

decreases  with  y  and  attains  its  minimal  quantity  ~  5 .  x  10  1  at  y  «  70  ^  4.4A  (blue  line). 

This  trend  is  due  to  the  fact  that  the  theoretical  weight  coefficient  (1.26)  is  the  ratio  of  scalar 
products,  which  are  integrals  over  the  half-infinite  interval  0  <  y  <  oo .  However,  in  calculations 
of  the  weight  coefficient  the  integration  is  perfonned  over  the  finite  interval  0  <  y  <  y  ;  i.e.,  the 

exponentially  decaying  tail  of  the  integrant  is  ignored.  This  leads  to  the  computational  error, 
which  decreases  as  the  upper  boundary  of  computational  domain  increases. 

Further  increase  of  y  leads  to  the  increase  of  computational  error  that  is  illustrated  in  Fig.  3  by 
the  magenta  line  for  y  =  90 .  The  minimal  error  depends  on  accuracy  of  the  mean-flow  profiles 
and  eigenfunctions,  which  are  used  in  calculations  of  the  scalar  products. 

To  verify  robustness  of  the  module  RECEPT  for  3-D  disturbances,  we  calculated  the  weight 
coefficient  for  the  disturbance  field  comprising: 

•  The  fast  acoustic  wave  of  the  angle  of  incidence  0(/  =  30°  and  various  sweep  angles 

0°  <  0_  <40°.  Its  phase  velocity  is  Ca  =  —  —  1  4 - - - ,  and  the  z- 

aa  Me  cos  Qy  cos  0Z 

component  of  wavenumber  is  f3a  =  aa  I  a  n  0  ; 

•  The  3-D  second-mode  wave  with  (3  =  0o  and  the  weight  coefficient  c2  =  1 
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The  disturbance  frequency  parameter  is  F  =  1 .3  x  10-4 ;  the  Reynolds  number 
R  =  ,JU*x*  /  v*e  =  10'! .  The  computational  grid  contains  501  equally  spaced  knots  with  the 
upper  boundary  y  =  50  . 

The  real  and  imaginary  parts  of  the  normal- wave  eigenvalue  a2 ((-),  )  are  shown  in  Figs.  4a, b. 
The  growth  rate  —  a  is  maximal  at  (-).  =  0  that  is  typical  for  the  second  mode.  Figure  5 
compares  the  actual  weight  coefficient  (black  line)  with  the  weigh  coefficient  c2  predicted  by 
the  RECEPT  (red  line).  The  computational  error  is  smaller  than  0.5%  at  0.  =  0  and  decreases 
with  the  sweep  angle. 


Fig.  la.  Mean-flow  velocity  profde;  Me  =  6  ,  Tw  =  T  ; . 
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T 


Fig.  lb.  Mean-flow  temperature  profile;  Mf  =  6 ,  Tw  =  T  ; . 


M  =6,  T  =7.05T 

e  w  e 


Fig.  2.  The  real  part  of  pressure  distributions  for  acoustic  wave  Fn(y),  the  second  mode  F 2(y) , 
and  the  disturbance  field  F n(y)  =  F a(y)  +  F 2{y) ;  R  =  103,  F  =  1.3  x  1CT4 . 
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Fig.  3.  The  weight  coefficients  predicted  by  RECEPT  at  various  upper  boundaries  y] . 


©z,  deg 

Fig.  4a.  Real  part  of  the  second-mode  eigenvalue  a  as  a  function  of  the  sweep  angle  0Z; 
M  =6,  T '  =T,  -R  =  103,  F  =  1.3  x  10~4 . 

e  w  ad 
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Fig.  5.  The  weight  coefficient  of  the  second  mode  predicted  by  the  RECEPT  for  various  sweep 
angles  Qz ;  Me  =  6 ,  Tw  =  Trf ,  R  =  103,  F  =  1.3  x  1CT4  . 
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5.  Receptivity  coefficients  for  2-D  second  mode  generated  by 
3-D  acoustic  waves  interacting  with  local  3-D  roughness 
element 

5.1  Problem  formulation  and  basic  relations 

Consider  the  boundary  layer  on  a  flat  plate  radiated  by  3-D  acoustic  wave  as  schematically 
shown  in  Fig.  6.  The  longitudinal  x* ,  normal  y* ,  transversal  C  coordinates  and  time  t*  are 

made  nondimensional  using  the  boundary-layer  length  scale  b  =  yjv*ex*  /  U*  and  U*  as 
(. x,y,z )  =  (x*,  y*,  z*) / 8* ,  t  =  t*U*  /  8* .  With  the  assumptions  that  the  wall  temperature  is 

constant  and  the  interaction  of  the  boundary  layer  with  the  inviscid  flow  is  weak,  the 
nondimensional  mean-flow  velocity  U ,  temperature  T ,  and  pressure  P  are  approximated  by 
the  compressible  Blasius  profiles 

U  =  U*  fU*  =  U{rj),  T  =  T* /T*  =  T(rj),  P  =  P*  /  P*  =  1,  (5.1) 

where  rj  =  y*  /  yj v*ex *  /  U*e  . 

At  the  distance  x*  from  the  leading  edge,  there  is  3-D  local  roughness  element  of  the  shape 

Vw\x  ,z)  =  h  o{x  -  x0)8{z  -  z0) ,  (5.2) 

where  8  is  Dirac  delta  function,  h*  is  roughness  height. 

A  three-dimensional  external  acoustic  wave  of  the  angular  frequency  uj  =  uf  8*  /  U*  and  wave- 
numbers  aa  =  a*a8* ,  j3a  =  /3*a8*  incidents  upon  the  plate  as  shown  in  Fig.  6.  The  nondimensional 
pressure  of  the  incident  wave  is  expressed  as 

P  —  £a  exp {iaax  +  i(3az  +  ikaincy  -  iut )  +  (c.c) .  (5.3) 

where  ea  is  acoustic-pressure  amplitude.  The  vertical  wave-number  component,  kamc,  is 
determined  from  the  acoustic  dispersion  relation  k  =  k  (a.  u.  M) ,  where 

ka  =  yj(M’  —  l)(aa  —  «,)(«„  —  aj  ,  Re(A:„)  >  0.  For  13 g  real  and  fixed,  the  branch  points  a12 
of  the  function  Kip)  correspond  to  slow  and  fast  acoustic  waves  with  the  wave  front  being 
perpendicular  to  the  wall;  i.e.,  the  angle  of  incidence  is  Q(/  =  tan_1(A:a  / aa)  =  0,  whereas  the 
inclination  angle  is  ©.  =  tan_1(/3a  / aa)  with  ka(a12)  —  0 .  The  wave  numbers  a  and  the 
correspondent  phase  speeds  C  =  u>  /  a  can  be  written  in  tenns  of  02 
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a 


1,2 


cuMcosO.  „  .  1 

_ -  /  '  -  I  x  _ 

M  cos  0.  =|=  1  ’  1,2  M  cos0  . 


(5.4) 


Hereafter  we  assume  that  M  cos0,  >  1 .  The  ranges  aa  >  a  and  aa  <  a  are  relevant  to  slow 
and  fast  acoustic  waves  with  non-zero  angle  0/; ;  the  wave-numbers  of  these  ranges  are 
expressed  as 


coM  cos0,  cos0, 

a0  = - ^  *  ,  (5.5) 

M  cos0;  cos©^  =F  1 

where  minus  (plus)  corresponds  to  slow  (fast)  acoustic  waves. 

The  incident  acoustic  wave  interacts  with  the  roughness  element  and  excites  boundary-layer 
modes  of  frequency  oo .  Hereafter  we  focus  on  excitation  of  2-D  second  mode,  which  is 
dominant  instability  at  sufficiently  high  Mach  numbers.  This  mode  is  given  by  the  vector- 
function 


'¥jx,y,z,t)  =  F0(x,y,z)exp{-iut)  +  (c.c). 


(5.6) 


Downstream  from  the  roughness  element,  x  >  x  ,  the  second-mode  amplitude  is  expressed  as 


x 

F2  =  heac(x0)A2(x,y)ex.p(iS2),  S2  =  Ja2(x,x)dx,  a2  =  ol2(x,u)  . 


(5.7) 


where  h  =  h*  /  6* ,  A2  is  eigenfunction  nonnalized  by  the  condition  H24(x,0)  =  1  (the  pressure 
amplitude  on  the  plate  surface  is  heac );  c  is  the  receptivity  coefficient,  which  depends  on  the 
roughness  element  locus  xQ  and  parameters  of  the  acoustic  wave.  Using  the  similarity  variables 
R.  7]  we  can  express  the  disturbance  amplitude  in  the  form 


F2  =  heac(RQ )A2(R,rj)  exp (iS2) ,  S2 


R 

2  J a2(R,u)dR ;  R>  RQ. 


As  shown  in  Section  1  the  receptivity  coefficient  is  expressed  as 


(5.8) 


c{R0)  =  qp(oiw0,(3w), 


(5.9) 
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where  q  is  given  by  Eq.  (1.14).  The  roughness-shape  Fourier  transform  p(aw,/3w)  is  given  by 
Eq.  (1.17).  For  the  roughness  shape  (5.1),  we  have 

+oo  +oo 

P(aw>Pw)  =  f  dz'  f  dx,6(x')6(z,)exp(-iawx'  -  i(3wzf)  =  1,  (5.10) 

— OO  —OO 

where  x  =x  —  x  ,  z  =  z  —  zQ  are  local  coordinates.  The  wave-numbers  ctw0,/3w  satisfy  the 
resonance  conditions  (1.23),  which  for  2-D  second  mode  with  (3  =  0  are  expressed  as 

awo  =  a2(xo /?„+/?„  =°.  (5.11) 

Substituting  (5.10)  into  (5.9)  we  obtain 

c(Ra)=q(R0,et,e„u).  (5.12) 

Using  the  receptivity  coefficient  (5.12)  we  calculate  the  initial  amplitude  of  the  second  mode  as 

V2(Ra,ri)  =  hES(R„)\(R0,r,) .  (5.13) 

In  particular,  the  initial  pressure  amplitude  on  the  plate  surface  is  given  by 

p2(it0,0)  =  FJR0, 0)  =  heac(R„).  (5.14) 


The  receptivity  coefficients  are  calculated  using  the  module  RECEPT  described  in  Section  2. 

5.2.  Numerical  results 

Consider  the  case:  M  =  6,  Prandtl  number  Pr=0.72,  specific  heat  ratio  7  =  1.4,  the  local 
temperature  T*  =  57.56  K,  the  wall  temperature  ratio  Tw  =  T*  / T*  =  5.212.  The  viscosity- 
temperature  dependency  is  given  by  Sutherland  formula 


Mr) 


(1  +  S)  z^-,3/2 

(T  +  S) 


(5.15) 


where  S  =  110./  T*  for  air  temperature  measured  in  Kelvin  and  fi  =  \i  j  /i*  is  nondimensional 
dynamic  viscosity.  The  ratio  of  the  bulk  viscosity  to  the  dynamic  viscosity  is  fiv/ n  =  0.8  (the 

corresponding  coefficient  in  the  stability  equations  is  e  =  —  —  =  1.2  ).  This  regime  was  used  in 

2  fi 


Ref.  [40]  for  parametric  studies  of  receptivity  to  2-D  acoustic  waves  interacting  with  2-D 
roughness  (waviness  or  hump). 


33 


ISTC  2393 


The  incident  acoustic  wave  is  specified  at  the  upper  boundary  r/,(/)  =  30^2  =  42.426,  which  is 

located  at  the  distance  from  the  wall  =2.5  of  the  boundary-layer  thickness  (see  Fig.  7a).  The 
disturbance  frequency  corresponds  to  the  frequency  parameter  F  =  8.x  10  5 ,  which  was 
considered  in  Ref.  [40]  for  2-D  case. 

The  mean  flow  profiles  U{rj)  and  T{rj)  are  shown  in  Figs.  7a  and  7b.  A  small  positive 
temperature  gradient  dT  /  dr/  on  the  wall  indicates  that  the  plate  surface  temperature  is  slightly 
less  than  the  adiabatic  wall  temperature. 

Distributions  of  Rea2(i?),  the  spatial  growth  rate  a  ~  —  1m  a  (R)  and  the  phase  speed 
C2  =  Re(cu  /  a  )  are  shown  in  Figs.  8a,b,c.  The  point  on  the  lower  neutral  branch  corresponds  to 
R  =  1000 .  The  non-monotonic  behavior  of  C^(R)  is  associated  with  the  presence  of  the 
spectrum  branchpoints  in  the  complex  R -plane  [21,23,40]. 

Distributions  of  the  receptivity  coefficient  modulus  |c(0J/,0z)|  for  the  roughness  element  loci: 
R0  —  500  (upstream  from  the  neutral  point),  R  =  1000  (at  the  neutral  point)  and  i?()  =  1500 

(downstream  of  the  neutral  point)  are  shown  in  Figs.  9a,b,c  for  fast  acoustic  waves.  The 
receptivity  coefficient  monotonically  decreases  with  the  inclination  angle  0 , ;  i.e.,  2-D  acoustic 
waves  interacting  with  3-D  local  roughness  element  provide  maximum  excitation  of  2-D  second 
mode.  On  the  contrary,  the  functions  |c(0y )|  at  0  =  fixed  have  maximums.  For  0.  =  0  (2-D 
acoustic  waves),  these  maximums  are: 

•  |c|  rj  7.  x  10~4  at  0^  rj  42°  and  RQ  =  500 

•  \c\  10.  x  10~4  at  0^  rj  38°  and  R(I  =  1000 

•  \c\  ~  6.  x  10~4  at  0^  Rj  48°  and  RQ  =  1500 

The  receptivity  coefficient  has  maximum  near  the  neutral  point  where  the  most  effective  angle  of 
incidence  is  close  to  its  minimum.  However,  the  receptivity  maximums  at  R  =  500  and  1500 

are  wider  than  that  at  R  =  1000  .  This  trend  may  compensate  the  maximum  in  the  neutral  point 
vicinity  if  the  incident  acoustic  waves  have  broad  distributions  with  respect  to  angles  0(/  and 

a- 

Numerical  results  for  slow  acoustic  waves  are  shown  in  Figs.  10a,b,c  for  the  same  loci  of 
roughness  element.  Similar  to  the  case  of  fast  acoustic  waves,  the  receptivity  coefficient 
monotonically  decreases  with  the  inclination  angle  0.  and  has  maximum  versus  the  angle  of 
incidence  0,, .  For  2-D  acoustic  waves  (0.  =  0  )  these  maximums  are: 

•  |c|  Ri  14.  x  10~4  at  rj  52°  and  RQ  =  500 

•  |c|  ~  16.  x  10~4  at  ©^  ~  52°  and  RQ  =  1000 
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•  |c|  «  22.  x  1(T4  at  <dy  «  30°  and  R{]  =  1500 

As  contrasted  to  the  case  of  fast  acoustic  waves,  the  receptivity  maximum  increases  downstream, 
while  its  width  is  non-mono  tonic:  it  is  essentially  larger  near  the  neutral  point  R  —  1000  .  Note 

also  that  slow  acoustic  waves  are  more  effective  than  fast  waves:  for  R  =  1000 , 
max  |c(@J/,0_, slow  waves)|  =1.6 max  fast  waves) | . 

Figures  lla,b  and  12a, b  show  that  the  real  and  imaginary  parts  of  the  receptivity  coefficients 
oscillate  versus  0 ;  i.e.,  the  initial  phase  of  the  second-mode  strongly  depends  on  the  angle  of 

incidence.  If  acoustic  waves  are  widely  distributed  with  respect  to  0  ,  these  oscillations  may 
lead  to  significant  reduction  of  the  second-mode  initial  amplitude. 


side  view 


top  view 


Fig.  6.  Schematics  of  incident  acoustic  wave  interacting  with  local  roughness  element. 
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Fig.  7a.  Profile  of  mean-flow  velocity. 
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Fig.  7b.  Profile  of  mean-flow  temperature. 
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Fig.  8c.  The  second-mode  phase  speed  C2  =  Re(cu  /  a  )  as  a  function  of  Reynolds  number  R. 
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Fig.  9a.  The  receptivity  coefficient  modulus  |c(09,0,)|  at  the  roughness  element  locus 
R  =  500  (upstream  from  the  neutral  point);  fast  acoustic  waves;  angles  are  given  in  degrees. 
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Modulus  of  c(0y,0z),  fast  acoustic  wave,  F=8.E-5,  R=1000 
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Fig.  9b.  The  receptivity  coefficient  modulus  |c(0,/?02)|  at  the  roughness  element  locus 
R0  =  1000  (neutral  point);  fast  acoustic  waves;  angles  are  given  in  degrees. 
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Fig.  9c.  The  receptivity  coefficient  modulus  c(0w,0.)  at  the  roughness  element  locus 
R0  =  1500  (downstream  from  the  neutral  point);  fast  acoustic  waves;  angles  are  given  in 
degrees. 
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Modulus  of  c(0y,©z),  slow  acoustic  waves,  F=8.E-5,  R=500 
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Fig.  10a.  The  receptivity  coefficient  modulus  |c(0y,0j|  at  the  roughness  element  locus 
Rq  =  500  (upstream  from  the  neutral  point);  slow  acoustic  waves;  angles  are  given  in  degrees. 
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Modulus  of  c(0y,0z),  slow  acoustic  waves,  F=8.E-5,  R=1000 
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Fig.  10b.  The  receptivity  coefficient  modulus  |c(0J/,0.)|  at  the  roughness  element  locus 
R  =  1000  (neutral  point);  slow  acoustic  waves;  angles  are  given  in  degrees. 
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Modulus  of  c(0  f0z),  slow  acoustic  waves,  F=8.E-5,  R=1500 
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Fig.  10c.  The  receptivity  coefficient  modulus  |c(09,0J|  at  the  roughness  element  locus 
R  =  1500  (downstream  from  the  neutral  point);  slow  acoustic  waves;  angles  are  given  in 
degrees. 


Fig.  11a.  The  real  part  of  c(0  ,0J  at  the  roughness  element  locus  R0  =  1000  (neutral  point); 
fast  acoustic  waves;  angles  are  given  in  degrees. 
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Fig.  lib.  The  imaginary  part  of  c(0  ,0  J  at  the  roughness  element  locus  R0  =  1000  (neutral 
point);  fast  acoustic  waves;  angles  are  given  in  degrees. 
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Fig.  12a.  The  real  part  of  c(0y,0j  at  the  roughness  element  locus  R0  =  1000  (neutral  point); 
slow  acoustic  waves;  angles  are  given  in  degrees. 
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Fig.  12b.  The  imaginary  part  of  c(©  ,0  J  at  the  roughness  element  locus  R0  =  1000  (neutral 
point);  slow  acoustic  waves;  angles  are  given  in  degrees. 
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6.  Receptivity  coefficients  for  3-D  Tollmien-Schlichting  waves 
generated  by  acoustic  waves  interacting  with  local  3-D 
roughness  element 

6.1.  Problem  formulation  and  basic  relations 

The  problem  is  formulated  similar  to  the  case  of  second-mode  excitation  discussed  in  Section  5. 
Consider  the  boundary  layer  on  a  flat  plate  radiated  by  3-D  acoustic  wave  as  schematically 
shown  in  Fig.  6.  The  nondimensional  mean-flow  velocity  U ,  temperature  T ,  and  pressure  P 
are  approximated  by  the  compressible  Blasius  profiles  (5.1).  At  the  distance  x*  from  the  plate 

leading  edge,  there  is  3-D  local  roughness  element  of  the  shape  (5.2).  The  nondimensional 
pressure  of  3-D  incident  acoustic  wave  is  given  by  (5.3).  This  wave  interacts  with  the  roughness- 
induced  disturbance  and  excites  3-D  Tollmien-Schlichting  (TS)  wave,  which  is  dominant 
instability  at  moderate  supersonic  Mach  numbers.  This  mode  is  given  by  the  vector- function 

xVTS{x,y,z,t)  =  FTS(Xly)exp(if3TSz  -  iut)  +  (c.c) .  (6.1) 

Downstream  from  the  roughness  element,  x  >  xQ,  the  TS-mode  amplitude  is  expressed  as 

X 

Fts  =  Kc(x0)Ara(z,y)exp(zSj,  STS  =  J  aTg(x,cv,/3)dx,  (6.2) 


where  h  =  h*  /  <*>*,  ATg  is  eigenfunction  normalized  by  the  condition  Aj,  (x,  0)  =  1 .  The 
receptivity  coefficient  c  depends  on  the  roughness  element  locus  xQ  and  parameters  of  the 
acoustic  wave.  Using  the  similarity  variables  f?,  r]  we  express  the  disturbance  amplitude  in  the 
fonn 


Fre  =  h£ac(R0)ATS(R^)^P{iSTS) ,  STS  =  2 J aTS(R,u,/3)dR;  R>  R0.  (6.3) 

Rq 

The  receptivity  coefficient  is  c(RQ)  =  qp(awQ,  (5W) ,  where  q  is  given  by  Eq.  (1.14).  The 
roughness-shape  Fourier  transform  p(aw ,  (3 w )  is  given  by  Eq.  (1.17).  The  wave-numbers  a w0,(3w 
satisfy  the  resonance  conditions  (1.23),  which  for  TS-wave  are  expressed  as 

«ro0  =  -  aa )  >  =  ±(/?ra  -  Pa)  ■  (6-4) 

Then  the  receptivity  coefficient  is 
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<R0) 


q+(Ra,Qy,Qz,u)  +  q_(R  ,Q  ,Qz,u) , 


(6.5) 


where  the  subscript  “±”corresponds  to  the  ±  component  of  (6.4).  Using  (6.5)  we  calculate  the 
initial  amplitude  of  the  TS  mode  as 

Tts(r o>v)  =  h£aC(R0)ATS(R0,v)  •  (6-6) 


In  particular,  the  initial  pressure  amplitude  on  the  plate  surface  is  given  by 

Pts(R  o’  °)  =  FtsSRv  °)  =  h£ac(R0 )  •  (6-7) 

The  receptivity  coefficients  are  calculated  using  the  module  RECEPT  described  in  Section  2. 

6.2  Numerical  results 

Consider  the  case  related  to  the  receptivity  experiments  of  Graziosi  and  Brown  [41]:  M  =  2.98  , 
Prandtl  number  Pr=0.72,  specific  heat  ratio  7  =  1.4 ,  the  stagnation  temperature  T*  =  278  K, 

the  plate  temperature  corresponds  to  the  adiabatic  wall,  T*  =  T*  .  The  viscosity-temperature 

dependency  is  given  by  Sutherland  formula  (5.15).  The  bulk  viscosity  is  assumed  to  be  zero, 
Hv  =  0 .  The  mean-flow  profiles  U(rj)  and  T[rf)  are  shown  in  Figs.  13a, b.  The  boundary-layer 
thickness  is  approximately  77  9.2  . 

The  incident  acoustic  wave  is  specified  at  the  upper  boundary  77  =  15>/2  21.2,  which  is 

located  at  the  distance  from  the  wall  =  2.3  of  the  boundary-layer  thickness  (see  Fig.  13a).  The 
disturbance  frequency  corresponds  to  the  frequency  parameter  F  =  1.4  x  10  " . 

Contour  plots  representing  the  distributions  of  ReaTS(R,F  )  and  the  spatial  growth  rate 
<jts  =  —  Im aTS(R,F)  are  shown  in  Figs.  14a, b,  where  F^=/3*u*e/U*  so  that  the 
nondimensional  transversal  wave-number  is  (3  =  F.R .  In  Fig.  14b,  the  boundary  of  colored 
region  represents  the  neutral  curve  in  the  R  —  F  plane.  At  R  =  900  the  maximum  growth  rate 
is  observed  at  F,  rs  0.68  x  10~4  that  corresponds  to  the  TS  wave-front  angle 

Fts  =  tan_1(/3/Re(aTO))  71°.  Hereafter,  the  roughness  element  is  located  at  the  station 

R  =  620 ,  which  is  close  to  the  neutral  curve  tip  (see  the  vertical  black  line  in  Fig.  14b). 

6.2.1  Fast  acoustic  waves 

Figure  15a  shows  the  receptivity  coefficient  c  as  a  function  of  7\  for  the  fast  acoustic  wave  of 

the  angles  Qy  =  45°  and  0.  =  0.  In  the  vicinity  of  F s  =  4’6  ><  10  5 ,  the  imaginary  part  and 

modulus  of  c  rapidly  drop.  A  nature  of  this  behavior  is  clarified  in  Fig.  15b,  where  the 
receptivity  coefficient  is  plotted  as  a  function  of  the  angle  ipw  between  the  z  -axis  and  the  wave 

front  of  the  resonant  component  of  roughness-induced  disturbance,  ipw  =  lau  '(0v;/aw)  (see 
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Fig.  16).  The  rapid  change  of  receptivity  coefficient  corresponds  to  (p*™™  =  70.4°  (see  magenta 
vertical  line  in  Fig.  15b)  at  which  the  wave-front  of  resonant  waviness  is  directed  along  the  Mach 
wave.  In  this  case,  the  freestream  velocity  component  nonnal  to  the  wave  front  corresponds  to 
the  transonic  Mach  number  Mn  =  1 .  For  <pw  <  (fT™ ,  the  resonant  waviness  generates 
supersonic  disturbances  outside  the  boundary  layer  since  Mn  >  1 .  For  tpw  >  iflans ,  these 
disturbances  are  subsonic.  Transition  from  supersonic  to  subsonic  waviness  is  accompanied  by 
the  singular  behavior  of  receptivity  coefficient. 

Similar  trends  are  observed  in  the  case  of  3-D  fast  acoustic  wave  with  0=0  and  0.  =  45° 

V  z 

(see  Figs.  17a, b).  For  this  case,  more  details  on  transition  through  the  transonic  regime  are 
shown  in  Figs.  18-20.  The  waviness-induced  disturbance  is  given  by  (1.9).  Outside  the  boundary 
layer,  the  vector  Fti  is  represented  by  the  sum  of  four  vectors 

K  =  J2 ^ 1  exP(Aj?/)  ’  Re(Aj)  <  o ,  y  -►  oo .  (6.8) 

j- 1 

One  of  these  vectors,  say  ,  is  relevant  to  inviscid  disturbances  of  acoustic  type;  the  others  are 
associated  with  vorticity  and  entropy  waves.  The  function  A ^(f^)  is  shown  in  Fig.  18.  For 

supersonic  waviness,  Re(AJ  is  almost  zero  and  Im(A1)>0;  i.e.,  the  disturbance  oscillates 

outside  the  boundary  and  forms  Mach  waves.  For  subsonic  waviness,  this  disturbance 
exponentially  decays.  Transition  from  supersonic  to  subsonic  regime  is  accompanied  by  rapid 
changes  of  A L(F  ) .  In  the  vicinity  of  the  transonic  point  F}  rs  5.5  x  1CT5 ,  where  X1  is  practically 

zero,  the  vector-function  F w(y)  varies  sharply  as  shown  in  Figs.  19,  20. 

Figures  15b  and  17b  show  that  the  modulus  of  receptivity  coefficient  in  the  case  of  supersonic 
resonant  waviness  is  significantly  larger  than  in  the  case  of  subsonic  waviness.  This  seems  to  be 
due  to  weak  penetration  of  subsonic  waviness-induced  disturbance  into  the  outer  inviscid  flow. 

Figure  21  shows  the  contour  plot  of  the  receptivity  coefficient  modulus  as  a  function  of  the 
acoustic -wave  angles  0,,  and  0  , .  The  transversal  wave-number  parameter  of  TS  wave  is  fixed, 

F  —  6.8  x  10  5 ,  and  corresponds  to  TS-wave  of  maximum  downstream  growth  (see  Fig.  14b). 
For  small  angles  of  incidence  0  <  5°,  the  receptivity  coefficient  has  maximum  versus  0. .  At 
Qy  =  0  this  maximum  is  close  to  the  inclination  angle  of  TS  wave  tpTS  rj  71° .  For  larger  angles 
of  incidence,  |c|  monotonically  increases  with  0Z.  For  02  >85°,  the  function  |c(0,;)|  has 
rnaxinium  |  c  9.  x  10  4  in  the  vicinity  of  0^  «  18° . 

6.2.2  Slow  acoustic  waves 

Figures  22a, b  show  the  functions  c(F ' )  and  c(<pw)  for  the  slow  acoustic  wave  of  0  =  45°  and 
0.  =  0.  The  receptivity  maximum  lies  in  the  range  of  subsonic  resonant  waviness.  Its  value, 
lc|max  ~  6.25  x  10-4 ,  is  close  to  the  receptivity  maximum,  |c|max  th  6.  x  10^4 ,  for  the  fast  acoustic 
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wave  of  the  same  angles  0y,02  (see  Figs.  15a, b).  Similar  to  the  case  of  fast  acoustic  waves,  the 
receptivity  coefficient  is  singular  at  ipw  =  70.4°  corresponding  to  the  transonic  resonant 
waviness. 

Figures  23a, b  show  the  functions  c(F  )  and  c(ipw)  for  the  slow  acoustic  wave  of  0  =  0  and 
0.  =  45° .  In  this  case,  there  are  tow  singular  points  associated  with  the  resonant  transonic 
waviness  of  <pw  =  4=70.4°  .  In  Fig.  24  a  contour  plot  of  the  receptivity  coefficient  modulus  is 
shown  as  a  function  of  0  and  0Z.  As  in  the  case  of  fast  acoustic  waves  (Fig.  21),  the 
transversal  wave-number  parameter  Ffj  =  6.8  x  10  5  corresponds  to  TS-wave  of  maximum 
downstream  growth.  The  red  dashed  lines  indicate  0  and  02  at  which  the  resonant  waviness  is 
transonic.  Between  these  lines  the  receptivity  coefficient  reaches  the  maximum 
Iclmax  ~  2.8  x  103  at  0y  fa  52°  and  0.  «  38°.  This  maximum  is  approximately  3  times  larger 
than  |c|max  for  fast  acoustic  waves  (see  Fig.  21). 
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Fig.  13a.  Profile  of  mean-flow  velocity,  boundary-layer  thickness  77  rj  9.2  ,  M  =  2.98  . 
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Fig.  13b.  Profile  of  mean-flow  temperature  on  adiabatic  plate,  M  =  2.98  . 
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Fig.  14a.  The  real  part  of  TS  eigenvalue  Re  aTS  as  a  function  of  Reynolds  number  R  and  F, . 
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Fig.  14b.  The  TS  growth  rate  aTS  =  —  Im  arg  as  a  function  of  Reynolds  number  R  and  F  . . 


Fig.  15a.  Receptivity  coefficient  as  a  function  of  F  for  fast  acoustic  wave  of  0  =  45°, 

P  y 

Oz  =  0;  R  =  620,  F  =  1.4  xlO"5. 
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Fig.  15b.  Receptivity  coefficient  as  a  function  of  <pw  for  fast  acoustic  wave  of  Q?/  =  45°, 
e2  =  0 ;  R  =  620 ,  F  =  1.4  x  1(T5 . 


Fig.  16.  Orientations  of  the  wave-front  of  resonant  waviness  with  respect  to  Mach  waves. 
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Fig.  17a.  Receptivity  coefficient  as  a  function  of  F  for  fast  acoustic  wave  of  0  =0, 

P  " 

0.  =  45° ;  R  —  620 ,  F  —  1.4  x  10~5 . 


9 

’w 


Fig.  17b.  Receptivity  coefficient  as  a  function  of  ipw  for  fast  acoustic  wave  of  0  =  0, 
0.  =  45° ;  R  —  620 ,  F  —  1.4  x  10^5 . 
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Fig.  18.  The  function  A 1(F  )  for  the  case  of  fast  acoustic  wave  of  Qy  =  0,  <dz  =  45° ,  R  —  620 , 
F  =  1.4  x  10~5 . 
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Fig.  19.  Distributions  of  the  x  -component  velocity  induced  by  the  wavy  wall  in  the  vicinity  of 
the  transonic  point  Fj  &  5.5  x  1CT5 . 
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Fig.  20.  Distributions  of  the  y  -component  velocity  induced  by  the  wavy  wall  in  the  vicinity  of 
the  transonic  point  F  «  5.5  x  1CT5 . 
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Fig.  21.  Modulus  of  receptivity  coefficient  as  a  function  of  the  fast  acoustic-wave  angles  0  and 
02 ;  R  =  620 ,  F  =  1.4  x  1(T5 ,  Fp  =  6.8  x  1(T5 . 
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Fig.  22a.  Receptivity  coefficient  as  a  function  of  F  for  slow  acoustic  wave  of  0  =  45° , 

p  y 

Qz  =  0 ;  R  =  620 ,  F  =  1.4  x  1(T5 . 


Fig.  22b.  Receptivity  coefficient  as  a  function  of  ipw  for  slow  acoustic  wave  of  0  =  45°, 
0  =  0 ;  R  =  620 ,  F  =  1.4  x  10~5 . 
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Fig.  23a.  Receptivity  coefficient  as  a  function  of  F  for  slow  acoustic  wave  of  0=0, 

p  * 

<c)z  =  45° ;  R  —  620 ,  F  —  1.4  x  10~5 . 


Fig.  23b.  Receptivity  coefficient  as  a  function  of  tpw  for  slow  acoustic  wave  of  0  =  0 , 
0S  =  45° ;  R  =  620 ,  F  =  1.4  x  10~5 . 
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Fig.  24.  Modulus  of  receptivity  coefficient  as  a  function  of  the  slow  acoustic-wave  angles  0?; 
and  Q;R  =  620 ,  F  =  1.4  x  10~5 ,  F  =  6.8  x  10~5 ;  dashed  red  lines  shows  0  and  0.  of 

p  y 

transonic  resonant  waviness. 
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7.  Recommendations  for  design  of  receptivity  experiments  in 
quiet  wind  tunnels 

Receptivity  problem  depends  on  many  parameters  including: 

•  Freestream  parameters 

•  Mean  flow  parameters 

•  Spectral  characteristics  of  acoustic  field  and  other  freestream  disturbances 

•  Shape  and  location  of  roughness  element 

To  reduce  the  number  of  independent  parameters  and  focus  on  basic  features  of  the  receptivity 
process,  experiments  should  be  conducted  on  bodies  of  simple  geometry,  namely,  on  a  flat  plate 
with  sharp  leading  edge  at  zero  angle  of  attack  (see  Fig.  6).  In  this  case:  the  bow  shock  is  weak 
and  its  effect  on  acoustic  field  is  minimal;  mean-flow  profiles  in  the  boundary  layer  are  easy  to 
simulate;  the  boundary-layer  disturbances  (first  and  second  modes)  are  well  known.  Even  for  this 
simple  configuration,  measurements  of  the  receptivity  coefficients  are  not  straightforward. 

The  major  challenge  is  associated  with  the  fact  that  it  is  difficult  to  separate  the  boundary-layer 
disturbances  induced  by  acoustic  waves  interacting  with  different  non-uniformities  of  the  mean 
flow.  For  example,  acoustic  waves  radiating  the  flat  plate  with  a  local  roughness  element 
generate  the  boundary-layer  modes  near  the  plate  leading  edge  and  in  the  roughness  vicinity.  If 
hot-wire  measurements  are  perfonned  in  the  boundary  layer  downstream  from  the  roughness 
element,  then  the  hot-wire  signal  includes  the  disturbances  from  the  both  regions  as  well  as  the 
acoustic  field  comprising  both  fast  and  slow  waves.  In  this  case,  it  is  difficult  to  identify  signals 
relevant  to  various  sources  of  receptivity. 

This  difficulty  may  be  resolved  using  acoustic  disturbances  of  relatively  short  duration  as 
sketched  in  Fig.  25.  If  the  acoustic  source  emanates  a  disturbance  comprising  few  wavelengths 
(wave  train),  then  the  acoustic  wave  train  and  the  wave  trains  of  boundary-layer  modes  are 
observed  in  different  time  intervals.  The  hot  wire  first  detects  the  train  of  fast  acoustic  waves, 
then  the  train  of  boundary-layer  waves  generated  by  the  fast  acoustic  waves  on  the  roughness 
element,  then  the  train  of  boundary-layer  waves  generated  by  the  fast  acoustic  waves  near  the 
plate  leading  edge.  Similar  sequence  of  signals  is  related  to  the  train  of  slow  acoustic  waves. 
Tracking  these  signals  it  is  feasible  to  measure  characteristics  of  every  wave  train  and  calculate 
the  correspondent  receptivity  coefficients. 

By  placing  the  acoustic  source  at  various  positions  with  respect  to  the  roughness  element  it  is 
possible  to  radiate  the  latter  by  acoustic  wave  trains  of  various  front  angles  0 ,  0.  (see  Fig.  26). 

This  makes  feasible  to  measure  the  receptivity  coefficients  as  functions  of  0I; ,  0.  and  verify  the 

correspondent  theoretical  dependencies. 

Our  theoretical  analysis  showed  that  receptivity  to  acoustic  waves  interacting  with  the  roughness 
element  strongly  depends  on  the  roughness  locus  and  shape.  For  experiments  on  excitation  of 
second-mode  disturbances,  which  are  dominant  in  high-speed  boundary  layer,  it  is  reasonable  to 
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focus  on  2-D  roughness  since  most  unstable  disturbances  are  two-dimensional.  This  may  be  a 
strip  mounted  on  the  plate  at  a  certain  streamwise  station  xs  as  schematically  shown  in  Fig.  27. 
The  theory  predicts  that  most  receptive  region  corresponds  to  the  lower  neutral  branch.  Placing 
the  strip  at  this  station  and  varying  the  free-stream  unit  Reynolds  number  one  could  observe  the 
receptivity  dependence  on  the  roughness  element  locus  relative  to  the  neutral  point.  Special 
attention  should  be  given  to  the  choice  of  the  strip  width  dx .  Our  theoretical  analysis  indicates 
that  receptivity  is  maximal  for  dx  corresponding  to  the  resonance  conditions;  i.e.,  the  strip  width 

should  be  chosen  as  dr  th  /  \an  ±  a „  ,  where  a-  a-  is  the  wave-number  of  second  mode  and 

x  /  |  Z  “  |  Z '  a 

acoustic  wave,  respectively.  For  more  accurate  estimates  of  the  strip  size  and  its  locus,  it  is 
recommended  to  conduct  calculations  of  receptivity  coefficients  for  particular  experimental 
conditions  using  the  module  RECEPT. 

For  boundary  layers  at  moderate  supersonic  speeds  (which  are  typified  by  the  experiments  of 
Graziosi  and  Brown  [41]  in  the  Mach=3  quiet  wind  tunnel  of  Princeton  University),  3-D  waves 
of  the  first  mode  (Tollmien-Schlichting  waves)  are  most  unstable.  Experiments  on  excitation  of 
TS  waves  should  be  focused  on  3-D  acoustic  waves  interacting  with  2-D/3-D  roughness 
elements.  It  is  reasonable  to  start  with  cylindrical  sources,  which  generate  predominately  2-D 
acoustic  waves,  and  the  periodic  3-D  roughness  structure  schematically  shown  in  Fig.  28.  The 
streamwise  size  dx  and  the  transversal  period  d.  of  this  structure  are  estimated  as 

dx  rj  2n  / \ax  ±  aa  | ,  2n  /  |/3X  | ,  where  (a] ,  (3^ )  are  wave-number  components  of  most 

unstable  TS  waves.  As  in  the  case  of  the  second-mode  excitation,  roughness  should  be  placed  at 
the  streamwise  station  relevant  to  the  lower  neutral  branch  of  the  first  mode. 

The  parametric  calculations  of  Section  6  show  that  the  receptivity  coefficient  is  singular  if  the 
roughness  is  transonic:  the  freestream  velocity  component  nonnal  to  the  front  of  roughness 
elements  corresponds  to  the  Mach  number  Mn  =  1  (see  Fig.  16).  It  is  interesting  to  verify  this 
theoretical  finding  by  testing  the  roughness  patterns  of  various  angularities  (pw  (from  supersonic 
to  subsonic). 

Since  receptivity  is  multi-parametric  problem  associated  with  delicate  physical  mechanisms, 
design  of  receptivity  experiments  should  be  guided  by  theoretical  modeling  of  experimental 
conditions.  This  may  be  done  using  the  receptivity  module  RECEPT.  With  such  a  combined 
theoretical  and  experimental  approach,  it  is  feasible  to  identify  and  investigate  various  features  of 
supersonic  boundary  layer  receptivity  by  performing  non-extensive  experimental  programs.  This 
helps  to  reduce  the  risk  and  cost  of  experiments,  especially  in  high-speed  wind  tunnels  where 
runs  are  quite  expensive. 
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Fig.  25.  Hot-wire  measurements  of  disturbances  generated  by  acoustic  wave  train  in  supersonic 
boundary  layer. 
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Fig.  26.  Acoustic  sources  at  various  positions  (red  and  black  points)  radiate  the  roughness 
element  by  wave  trains  of  various  front  angles  0?/ . 


free  stream 


side  view 


Fig.  27.  Receptivity  to  acoustic  waves  interacting  with  2-D  roughness  element  (strip). 
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Fig.  28.  Receptivity  to  acoustic  waves  interacting  with  3-D  periodic  roughness  element. 
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8.  Conclusions 

1.  Analysis  of  high-speed  boundary  layer  receptivity  to  3-D  external  waves  interacting  with  3-D 
roughness  element  or  waviness  was  perfonned  and  summarized  in  the  AIAA  Paper  2003- 
3137  (33rd  Fluid  Dynamics  Conf.  and  Exhibit,  23-26  June  2003,  Orlando  Florida) 

2.  The  FORTRAN  computational  module  RECEPT.FOR  was  developed  and  delivered  to  the 
end  user  Dr.  Roger  Kimmel  with  the  usage  instructions.  This  module  allows  for  calculations 
of  the  receptivity  coefficients  and  the  weight  coefficients  for  modes  of  discrete  spectrum  (TS 
waves,  second  mode  etc.). 

3.  It  was  shown  that  the  RECEPT  (with  the  option  providing  calculations  of  the  weight 
coefficient  for  a  particular  discrete  mode)  allows  for  extraction  of  the  nonnal  wave  of 
relatively  small  amplitude  (-0.1%)  from  the  disturbance  field.  With  this  tool,  the  user  can 
identify  nonnal  waves,  which  latent  in  a  “messy”  disturbance  field  at  early  stages  of 
instability.  This  is  of  particular  importance  for  analyses  of  DNS  results  and  experimental  data 
in  connection  with  receptivity  and  stability  problems. 

4.  The  module  RECEPT  was  applied  to  the  problem  of  high-speed  boundary-layer  receptivity  to 

acoustic  disturbances  interacting  with  3-D  roughness  element.  Namely,  the  receptivity 
coefficients  were  calculated  for  2-D  second-mode  excited  by  3-D  incident  acoustic  wave 
interacting  with  3-D  roughness  element  on  a  flat  plate.  These  calculations  were  conducted  for 
the  free-stream  Mach  number  6  at  flow  and  wall  temperatures  typical  for  quiet  wind  tunnels. 
It  was  found  that  the  receptivity  coefficient  monotonically  decreases  with  the  acoustic-wave 
inclination  angle  (~)z  for  both  fast  and  slow  incident  waves.  2-D  acoustic  waves  interacting 
with  3-D  local  roughness  element  provide  maximum  excitation  of  2-D  second  mode.  At 
0  .  =fixed,  the  modulus  of  receptivity  coefficient  has  maximum  versus  the  angle  of  incidence 
0  .  In  the  case  of  0.  =  0  (2-D  acoustic  waves)  this  maximum  is  observed  at  0^  40°  for 

fast  waves,  and  0(;  rs  50°  for  slow  waves.  Receptivity  to  slow  acoustic  waves  is  higher  than 

to  fast  waves.  The  initial  phase  of  the  excited  second  mode  strongly  depends  on  the  angle  of 
incidence.  This  may  significantly  reduce  the  second-mode  initial  amplitude,  if  acoustic 
waves  are  distributed  over  a  wide  range  of  0  . 

5.  The  receptivity  coefficients  were  calculated  for  3-D  Tollmien-Schlichting  waves  (first  mode) 
generated  by  acoustic  waves  interacting  with  local  roughness  element.  These  calculations 
were  performed  for  the  boundary  layer  on  a  flat  plate  at  Mach=3  free  stream  with  parameters 
relevant  to  the  receptivity  experiments  of  Graziosi  and  Brown  [41].  It  was  found  that  the 
receptivity  coefficient  is  singular  if  the  resonant  wall  waviness  is  transonic;  i.e.,  the  free- 
stream  velocity  component  nonnal  to  the  waviness  front  equals  the  speed  of  sound. 
Numerical  examples  show  that  maximum  of  the  receptivity  coefficient  associated  with  slow 
acoustic  waves  is  essentially  larger  than  that  associated  with  fast  waves. 

6.  Recommendations  for  design  of  receptivity  experiments  in  supersonic  quiet  wind  tunnels 
were  fonnulated.  It  was  recommended  to  conduct  receptivity  experiments  on  a  flat  plate  in 
order  to  focus  on  basic  features  of  receptivity  process  and  minimize  the  number  of 
independent  parameters.  It  was  pointed  out  that  the  major  challenge  of  receptivity 
measurements  is  associated  with  the  fact  that  it  is  difficult  to  sort  out  the  boundary-layer 
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disturbances  induced  by  acoustic  waves  interacting  with  different  non-uniformities  of  the 
mean  flow.  This  difficulty  may  be  resolved  using  acoustic  wave  trains. 

7.  Since  receptivity  is  multi-parametric  problem  associated  with  delicate  physical  mechanisms, 
design  of  receptivity  experiments  should  be  guided  by  theoretical  modeling  of  experimental 
conditions,  which  could  be  performed  using  the  module  RECEPT.  With  such  a  combined 
theoretical  and  experimental  approach,  it  is  feasible  to  identity  and  investigate  various 
features  of  supersonic  boundary-layer  receptivity  by  conducting  non-extensive  experimental 
programs.  This  helps  to  reduce  the  risk  and  cost  of  experiments,  especially  in  high-speed 
wind  tunnels  where  runs  are  quite  expensive. 


9.  Future  effort 

Another  mechanism  of  receptivity  is  associated  with  nonlinear  interactions  of  freestream 
disturbances.  Choudhari  and  Streett  [22]  pointed  out  that  the  bilinear  interaction  between  a  pair 
of  freestream  waves  could  be  important  in  excitation  of  unstable  disturbances.  Wu  [42] 
elaborated  this  mechanism  considering  suitable  convecting  gusts,  which  interact  with  sound 
waves  producing  a  forcing  that  has  the  same  time  and  length  scales  as  those  of  Tollmien- 
Schlichting  waves  in  the  neutral  point  vicinity.  The  sound-convective  gust  interaction  may  be 
dominant  in  boundary  layers  at  subsonic  and  low  supersonic  speeds. 

For  high-speed  boundary  layers  it  is  feasible  to  identify  pare  of  acoustic  waves,  which  satisfies 
the  resonance  condition  un  =  u>  ±  c ca2  and  the  synchronization  conditions  an  =  a  ±  aa2 , 

[3n  =  (3  ±  l3a2 .  Here  u>a,  aa,  13 a  -  frequency  and  wavenumber  components  of  acoustic  waves; 
un ,  an,  (3n  -  correspondent  characteristics  of  the  boundary-layer  mode  (this  could  be  the  first 
n  —  1  or  second  n  —  2  mode).  Under  these  conditions,  the  scale  conversion  is  achieved  without 
resorting  of  a  local  non-homogeneity  of  the  mean  flow  (roughness  element,  waviness  etc.),  that 
may  give  large  receptivity  coefficients  on  smooth  surfaces.  Moreover,  one  of  interacting  acoustic 
waves  may  belong  to  the  low-frequency  band,  which  corresponds  to  intense  disturbances 
especially  in  the  wind-tunnel  free  stream. 

Our  preliminary  estimates  show  that  the  bilinear  sound-sound  interaction  induces  the  boundary- 
layer  modes  of  the  amplitude  An0  ~  /l;  /  A1/ 1 ,  where  R  is  displacement  thickness  Reynolds 

number.  For  high-speed  boundary  layers,  the  Reynolds  number  is  R  ,  ~  104 .  Assuming  that  the 

8 

transitional  amplitude  A  tr  ~  10  and  the  amplification  factor  N  ~  6 ,  we  estimate  the  initial 
amplitude  as  An0  ~  10-2e-6  ~  10-5,  which  corresponds  to  the  acoustic-wave  amplitude 

Aa  ~  . /i'1  ~  VlO  =  0.1%.  Such  levels  of  acoustic  disturbances  may  be  generated  by 

turbulent  boundary  layers  on  forebody  surfaces  (in  flight)  or  nozzle  walls  (in  wind  tunnels).  This 
leads  to  expectation  that  the  bilinear  sound-sound  interaction  may  dominate  the  receptivity 
process. 
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To  our  knowledge  the  discussed  above  mechanism  has  not  been  studied  yet.  We  suggest  to 
perform  theoretical  studies  of  the  bilinear  sound-sound  interaction  and  incorporate  this 
receptivity  mechanism  into  the  computational  module  RECEPT. 
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Abstract 

Theoretical  study  of  hypersonic  boundary-layer  receptivity  to  acoustic  waves  interacting  with  wall  waviness 
(distributed  roughness)  or  a  hump  (local  roughness)  is  performed  using  a  combination  of  asymptotic  and 
numerical  techniques.  It  is  shown  that  strong  excitation  of  unstable  disturbances  occurs  in  local  regions  where 
the  forcing,  which  is  due  to  bilinear  interaction  of  external  waves  with  wall-induced  disturbances,  is  in 
resonance  with  the  boundary-layer  modes.  The  receptivity  coefficients  are  presented  for  two-dimensional 
second  mode  generated  by  acoustic  disturbances  of  various  frequencies  and  incident  angles  interacting  with 
distributed  and  local  roughness  on  a  flat  plate  at  Mach  6.  The  theoretical  model  can  be  used  to  predict  initial 
amplitudes  of  the  second  mode,  which  is  dominant  instability  in  essentially  two-dimensional  boundary  layers  at 
hypersonic  speeds. 


At 


C 

F  =  co* v*e /U*e2 


Nomenclature 

=  eigenvector(direct  problem) 
=  eigenvector(adjoint  problem) 
=  phase  speed 
=  frequency  parameter 


/ 

G 

g 

L 

M 


=mass-flow  disturbance 
=  receptivity  function 
=  roughness  shape-function 
=  streamwise  scale 
=  Mach  number 


*  *2 

p  =  p*  /( peU  ~ )  =  pressure  disturbance 
p  =  p*  /  P*  =  pressure  disturbance 
q  =  receptivity  factor 


t 


=  time 


T 

u,v 

n,v,w 

x,y,z 


=  mean  flow  temperature 
=  mean  flow  velocity 
=  velocity  disturbance 
=  Cartesian  coordinates 


x1  =  x"  /  Z*  =  slow  variable 

a,  P  =  wave-vector  components 


8  =  5"  IL' 

A* 


v 

CO 


=  small  parameter 
=  boundary-layer  thickness 
=  kinematic  viscosity 
=  angular  frequency 


0 

=  angle  of  incidence 

9 

=  temperature  disturbance 

P(*,y) 

=  amplitude  vector-function 

Subscripts 

0 

=  resonance  point 

a 

=  acoustic 

ad 

=  adiabatic 

b 

=  branch  point 

dis 

=  distributed 

e 

=  upper  boundary-layer  edge 

loc 

=  local 

=  wall 

Superscripts 

bar 

=  complex  conjugate 

asterisk 

=  dimensional 

1.  Introduction 

Important  mechanism  of  boundary-layer 
receptivity  is  the  interaction  of  acoustic  disturbances 
with  steady  disturbances  induced  by  spatial  variations 
of  boundary  conditions  on  the  wall  surface  such  as 
roughness/waviness,  local  heating/cooling  and 
suction/blowing.  Early  experimental* 1  and  theoretical2' 
4  studies  of  this  mechanism  were  performed  in  Russia 
for  subsonic  flows.  Their  results  were  summarized  in 
Ref.  [5].  An  asymptotic  theory  of  the  Tollmien- 
Schlichting  wave  excitation  due  to  interaction  of 
acoustic  waves  with  surface  irregularities  was 
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developed  in  Refs.  [6-8]  and  reviewed  in  Refs. 
[8-10].  Detailed  parametrical  studies  were 
represented  in  Refs.  [11-14],  The  theoretical 
predictions  are  in  good  agreement  with  the 
experiment15  conducted  on  a  flat  plate  in  low- 
speed  wind  tunnel. 

Receptivity  of  high-speed  boundary  layers  (in 
the  Mach  number  range  from  1.15  to  4.5)  to 
freestream  disturbances  interacting  with  surface 
irregularities  was  analyzed  by  Choudhari  and 
Streett16.  They  reported  that,  although  the 
Mach=4.5  boundary  layer  is  nearly  6  times 
thicker  than  in  the  Mach=0  case,  the  instability 
amplitudes  generated  by  a  wall  hump  of  fixed 
height  are  only  marginally  smaller  than  those  at 
low  speeds.  This  example  indicates  that 
receptivity  of  hypersonic  boundary  layer  to 
acoustic  plus  roughness  disturbances  may  play 
important  role  in  the  initial  phase  of  transition. 

In  this  paper,  we  discuss  theoretical  analysis 
of  hypersonic  boundary-layer  receptivity  to 
acoustic  waves  interacting  with  surface  waviness 
(distributed  excitation)  or  hump  (local  excitation) 
on  a  flat  plate.  In  Section  2,  we  formulate  the 
problem  and  show  that  its  main  features  are 
similar  to  the  case  of  receptivity  to  wall 
disturbances.17  In  Section  3,  we  discuss 

numerical  results  related  to  excitation  of  the 
second  mode  by  fast  and  slow  acoustic  waves  of 
various  angles  of  incidence.  In  Section  4,  we 
conclude  the  paper. 

2.  Problem  formulation  and  analysis 

Hereafter  we  consider  excitation  of  the 
boundary-layer  modes  by  acoustic  waves 

radiating  a  flat  plate  with  a  wavy  surface 
(distributed  receptivity)  or  a  hump  (local 

receptivity).  Analysis  is  conducted  using  the 
method17  for  receptivity  to  wall-induced 

disturbances. 

2.1  Distributed  receptivity  on  wavy  surface 

Consider  two-dimensional  laminar  boundary 
layer  of  perfect  gas  on  a  flat  plate  in  hypersonic 
freestream.  The  global  streamwise  length  L *  is 
assumed  to  be  much  larger  than  the  boundary 

layer  thickness  scale  8"  =  *Jv*eL*  IU*e  ;  i.e.  the 
ratio  s  =  S*  /  L*  is  small.  The  longitudinal  x* , 
normal  y * ,  transversal  z*  coordinates  and  time 
f  are  made  nondimensional  as  (x,y,z)  = 
(x*,  y*,3  *)/ 8" ,  t  =  fU'e  /  8" .  Introducing  the 
slow  variable  xx  =  x*  /  V  =  sx  we  can  specify 


the  nondimensional  mean-flow  velocity  (U,  V) ,  and 
temperature  T  as 

U  =  ur  /U :  =  U (Xj,y) ,  V  =  r  /u;  =  £V0(x„y)  , 

T  =  r  /T*  =  T(xx,y) ,  (1) 

The  wall  waviness  is  represented  in  the  form 
yjx,  y,  z)  =  h  exp (iawx  +  ij3wz)  +  ( c.c ) ,  (2) 


where  h  =  h"  IS"  is  nondimensional  amplitude, 
(aw,j3w)  is  wavenumber  vector.  It  is  assumed  that 
the  wavelength  is  of  the  order  of  the  boundary-layer 
thickness  A*  ( U"  =  0.991/*  at  y*  =  A* ).  With  h * 
being  smaller  than  the  viscous  sublayer  thickness,  the 
wavy  wall  generates  small  steady  disturbances 
governed  by  the  linearized  Navier-Stokes  equations.18 

An  external  acoustic  wave  of  the  angular 
frequency  co  and  the  wavenumber  ( arl,j3a ), 
incidents  upon  the  plate  as  schematically  shown  in 
Fig.  1.  Its  physical  quantities  are  expressed  in  a 
traveling-wave  form.  For  example,  the  pressure  is 

P  =  P a, inc  exp [i(aax  +  paz  +  ka  incy  -  wt)]  +  (c.c) .  (3) 


The  vertical  wave-number  component,  kainc,  is 
determined  from  the  acoustic  dispersion  relation 
kaJ„c=kai„c(aa,  a,  M),  where  M  is  freestream 

Mach  number.  The  incident  wave  is  refracted  by  the 
boundary  layer  and  reflected  by  the  wavy  wall. 
Interacting  with  the  waviness-induced  disturbance  the 
acoustic  wave  generates  boundary-layer  normal 
waves  including  unstable  modes.  Our  objective  is  to 
determine  the  receptivity  function  coupling 
characteristics  of  the  incident  acoustic  wave  and  the 
wall  waviness  with  the  normal  waves. 

Using  the  approach17  we  represent  the  unsteady 
flow  field  by  the  vector  function 


/  du  _  ^  dO  5w.r 
™(x,y,t)  =  (u,  —,v,p,  0,  —,w,—)  . 

dy  dy  dy 


(4) 


which  is  expressed  in  the  form 

'V  =  Q(Xj ,  y)  +  hFw  (x, ,  y)  exp  (iawx  +  ij3„,z)  + 

+£aFa  (Xj ,  y)  exp  (iaax  +  ij3az-  icot) 

+F(x,  y,  z)exp(-icot)  +  (c.c)  ^ 

The  first  term  of  (5)  describes  the  mean  flow  on  a 
smooth  surface;  the  second  term  corresponds  to  the 
stationary  disturbance  induced  by  the  wall  waviness; 
the  third  term  corresponds  to  the  acoustic  disturbance. 
Definition  of  the  amplitude  sa  depends  on 
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normalization  of  the  incident  acoustic  wave. 
Hereafter  the  vector-function  Fa(j)  is 
normalized  as  (pressure  amplitude  of  the  incident 
acoustic  wave)=l  at  the  upper  boundary-layer 
edge,  i.e.  ea  =  p*Jnc/P*  ■  The  fourth  term  of  (5) 

represent  disturbances  generated  in  the  boundary 
layer  due  to  nonlinear  interaction  of  the  acoustic 
wave  with  the  waviness-induced  perturbation. 

Substituting  (5)  into  Navier-Stokes  equations 
we  obtain  the  following  problems  written  in  a 
matrix-operator  form. 

For  disturbances  from  wavy  wall: 

H(y,dy,xx,to  =  0,aw,j3JFw  =  0  ,  (6a) 

Fwl =  Fw 3 (x, , 0)  =  0 , 

Fw5 (xi,0)  =  -^(x1,0),  Fw1  (x,,0)=0,  (6b) 

|FW|  <  oo  ,  y  — »  oo  ,  (6c) 

where  the  boundary  conditions  (6b)  are  obtained 
using  the  Taylor  series  expansion  of  the  no-slip 
conditions,  u  =  v  =  w  =  0 ,  and  the  wall- 
temperature  condition,  9  =  0  ,  at  y  =  yw. 

For  acoustic  disturbances: 

H(y,dy,xi,(o,aa,Pc) Fu  =  0  ,  (7a) 

Fa,=Fai=FaS=Fal=0,y  =  0,  (7b) 

(pressure  amplitude  of  incident  wave)=l,  y  — >  oo  . 


modes  Ft  (x, ,  y)  with  the  eigenvalues  ak  (x, ,  P,  to) 

F  =  hea  Yj  [Fao  (x, ,  y)  +  eFkl  (x, ,  y)  +  ■  ■  •]  x 

k 

exp (i£~lSk +ij3z) ,  (9a) 

Sk  =  J  ak  (x, ,  co,  P)dxi  .  (9b) 

If  the  nonlinear  term  of  (8a)  is  not  in  resonance 
with  the  boundary-layer  modes, 

ak  (x, ,  to,  P)  ^  aw  +  aa  and/or  P  #  Pw  +  Pa ,  then 
receptivity  is  exponentially  weak.  In  this  case,  we 
obtain  the  eigenvalue  problem  typical  for  linear 
stability  analyses 


-f-//oK=o, 

) 

(10a) 

(F*o)l,3,5,7  =  y  ~  0  , 

(10b) 

o,  y  -> 00  • 

(10c) 

Solutions  of  (10)  describe  the 

boundary-layer 

normal  waves  far  from  resonance  regions.  They  can 

be  expressed  as 

F,0  =ck(xl)Ak(x],y,ak), 

(11) 

where  Ak  are  eigenfunctions  normalized  by  a  certain 
condition  such  as  Ak4(x],0,ak)  =  1  .  The  amplitude 
functions  ^.(Xj)  are  unknown  on  this  stage. 

If  the  nonlinear  term  is  in  resonance  with  one  of 
the  normal  waves  (say  with  a  =  an )  at  the  point 


For  the  boundary-layer  disturbance: 

H(y,  8y ,  xx ,  sdH ,  z,  d: ,  oj)F  = 
hsfi{y, x, , Fb, ; <o, aw , P„ ,aa,P„ )Fa  x 
exp  [i(aw  +aa)x  +  i(Pw  +  Pa  )z] ,  (8a) 

Ft  =  F3  =  F5  =  F7  =  0 ,  y  =  0 ,  (8b) 

|F|  <  oo  ,  y  — >  oo .  (8c) 

The  right-hand  side  of  (8a)  describes  bilinear 
interaction  of  the  wall-induced  disturbance  with 
the  acoustic  wave;  the  higher  order  terms  ~  s2Ji , 
sclh2  are  neglected.  In  (6a),  (7a)  and  (8a),  the 
operator  FI  can  be  expressed  in  the  standard 

form  H  =  -^—-H0(y,xl',R,to,a,P)  ,  where  the 

dy 

8x8  matrix  H0  depends  on  the  mean-flow 
profiles  (1)  and  disturbance  parameters.19 

We  consider  a  partial  solution  of  the  problem 
(8a)-(8c)  as  a  decomposition  of  the  discrete 


Xi  x10  , 

a„ (xUi,oj,P)  =  aw  +  a.  and  P=Pw  +  Pa,  (12) 

then  the  acoustic  wave  effectively  generates  this 
normal  wave  due  to  the  sound-waviness  interaction 
near  the  resonance  point  x10 . 

We  consider  the  case  when  the  normal  wave  is 
unstable  and  the  resonance  point  is  close  to  the  lower 
neutral  branch.  Analyzing  the  governing  terms  of  (8a) 
we  conclude  that  the  resonance  excitation  occurs  in 
the  relatively  narrow  region  |xj  -x10|  =  0(s112)  with 

the  inner  variable  £  =  £~1/2(x, -x10)  =  0(1) .  We 
assume  that 

P  =  Pw  +  Pa  and  aw  +  aa  =  a0  +  ye 1/2 ,  (13) 

where  a0  =  an(xl0,to,P) ;  the  real  parameter  y 
characterizes  resonance  detuning.  Then,  the  inner 
solution  of  the  problem  (8a)-(8c)  is  expressed  as 
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F  =  sah[e~U2c(^)  A0  (y)  +  A,  (£  y) . . .]  x 
exp  (i£~'/2a0%  +  ipz) 


where  A0(y)  =  An(xl0,y,a0)  is  the  normal- 
wave  eigenfunction  at  the  resonance  point. 
Expanding  the  matrix  H0  in  the  vicinity  of 
x  =  x10  we  obtain  the  problem  for  A, 

~z  Ho(xio’ao) 

V<?y 

f^o. (x|(),«0) A(l  +  GFa  exp(iyp) ,  (15a) 

oxi 

(4W=o,  y  =  o,  (15b) 

|A,|->0,  J^00-  (15c) 

This  problem  is  similar  to  (2.18)  of  Ref.  [17]. 
It  has  non-trivial  solution  if  its  right-hand  side  is 
orthogonal  to  the  eigenfunction 
,y,a0)  of  the  adjoint  problem.  This 
condition  leads  to  the  equation  for  c(c) 


.  dc  dH0 

Ai  (*i<»«o)Ao  + 

dq  da 


This  amplitude  is  proportional  to 

(Z*  IS*)'12  =  s~'12 ,  whereas  the  resonant  region 

length  is  r=s~'l28*.  In  the  parallel  flow 
limit,  s  — »  0 ,  the  resonance  region  is  infinitely  large 
and  the  normal-wave  amplitude  tends  to  infinity.  In 
this  case,  the  upper  limit  of  the  disturbance  amplitude 
is  determined  from  analysis  accounting  for  nonlinear 
terms  of  higher-order.  Due  to  nonparallel  effects  the 
receptivity  mechanism  is  local  in  units  of  Z *  even  for 
distributed  roughness.  However,  it  can  be  treated  as  a 
distributed  mechanism  in  units  of  8 * .  Our  choice  of 
the  term  “distributed  receptivity”  is  a  matter  of 
convention. 

2.2  Local  receptivity  on  hump 

We  consider  a  span  wise  row  of  local  humps  on  a 
flat  plate  surface.  For  brevity  we  analyze  one 
spanwise  harmonic  of  this  periodic  structure,  and 
specify  its  shape  as 

yjx,  y,  z)  =  hg(x’)exp(ifiwz)  +  (c.c)  ,  (20) 

X  =  X  -  X0  =  s  '  (X[  -  x10 ) . 


uu 

~d^~ 

-ib^c  =  q  exp(iy^)  , 

(16a) 

b  = 

^L(r  ) 

,  vAio/> 

ax  j 

(16b) 

q  = 

i(B0,G(x1o)F<J) 

(16c) 

(B°’  da 

where  the  scalar  product  is  defined  as 

00  8  __ 

(B’AW  iuBjAjdy,  (17) 

0  7=1 

the  overbar  denotes  complex  conjugate  values.  If 
the  normal  wave  has  zero  amplitude  upstream 
from  the  resonance  region,  then  the  boundary 
condition  for  (16a)  is 

c(4)->0,  (18) 


The  problem  (16a)-(18)  is  similar  to  the 
problem  (2. 1 9)-(2.2 1 )  of  Ref.  [17].  Using  this 
similarity  we  obtain  the  normal-wave  amplitude 
downstream  from  the  resonance  region 


t-'  7  -1/2  2777 

F„  =  sjis  -q  —  exp 
V  b 


b  jy2  ^ 
K~2 ~bj 


exp 


■*i 

is  '  j  a„  (x,  )dx] 


,  Xj  >  XjQ  . 


A„(w)x 

(19) 


The  shape-function  is  represented  by  the  Fourier 
integral 

+oo 

g(x')  =  —  f  P(aj  exp (iawx')daw ,  (21) 

which  can  be  expressed  in  the  form 

1/2  -Ko 

g  =  —exp (iaw0x')  f  p(aw0  +  s"2y)e'r(dy  .  (22) 

2n  i 

Here  aw0  =  an(x10,a>,fl)  -aa  is  the  wave-number 
component  determined  from  the  resonance  condition 
(12);  y  characterizes  resonance  detuning  in 
accordance  with  (13).  Using  (19)  for  the  normal-wave 
amplitude  excited  by  one  harmonic  and  integrating 
over  y ,  we  obtain 


F„ 


=  sahqp{aM)\n  exp 


(  x'  > 

an(x')dx'  +  i(lz 

Vo  ) 


(23) 


Note  that  the  function  q(xl0,a0,aa,  fia)  does  not 
depend  on  the  hump  shape  g(x') .  If  this  function  is 
calculated  once  and  for  all,  then  we  can  evaluate  the 
normal-wave  amplitude  for  a  wall  hump  of  arbitrary 
shape  (within  the  limit  of  basic  assumptions). 

Equation  (23)  indicates  that  acoustic  waves 
interacting  with  a  local  hump  generate  normal  waves 
of  a  finite  amplitude  even  in  the  parallel  mean  flow. 
This  is  due  to  the  fact  that  a  small  portion  of  the 
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hump  spectrum,  (aw  -aw0)  =  0(s'/2) ,  is 
involved  into  the  resonance.  As  e  — >  0 ,  the 
resonant  excitation  increases  as  s~"2 ,  whereas 
the  active  spectral  band  is  narrowed  down  -  s'1 2 . 
Because  both  trends  compensate  each  other,  the 
normal-wave  amplitude  remains  finite. 

The  relations  (19)  and  (23)  hold  for  any 
combinations  of  freestream  disturbances  and 
wall  perturbations.  If  we  replace  the  acoustic 
wave  by  vortical  (entropy)  wave,  then  these 
expressions  provide  the  normal-wave  amplitude 
generated  by  vortical  (entropy)  wave  interacting 
with  the  surface  irregularity.  The  normal  wave 
can  also  be  of  different  nature:  Tollmien- 
Schlichting  wave,  second  mode,  cross-flow 
vortices,  Gortler  vortices  etc.  All  these  cases  can 
be  treated  using  the  following  approach: 

1 .  Solve  the  problem  (6a)-(6c)  for  wall-induced 
disturbance 

2.  Solve  the  problem  (7a)-(7c)  for  external 
(freestream)  disturbance 

3.  Solve  the  eigenvalue  problem  (10a)-(10c) 
and  the  adjoint  problem  to  obtain  the 
eigenfunctions  A „(x1,y)  and  BJ1(x1,>’) 

4.  Evaluate  the  scalar  products  of  (16c) 

5.  Predict  the  normal-wave  amplitude  using 
(19)  for  distributed  receptivity  or  (23)  for 
local  receptivity 

This  algorithm  can  be  incorporated  into  a 
computational  module,  which  may  be  linked 
with  a  standard  linear  stability  solver. 

2.3  Receptivity  near  spectrum  branch  points 

Foregoing  analysis  is  valid  for  “simple” 
discrete  spectrums;  i.e.  the  eigenvalues  of 
boundary-layer  modes  are  assumed  to  be 
essentially  different.  This  is  typical  for  boundary 
layers  at  subsonic  and  moderate  supersonic 
speeds  as  well  as  hypersonic  boundary  layers  on 
adiabatic  walls.  Flowever,  in  the  case  of  low  wall 
temperature  ratios,  the  first  and  second  modes 
have  the  branch  points,  which  are  close  to  the 
real  parameter  space.17,21'22  The  upstream  branch 
point,  x,  =  xlb ,  is  located  near  the  lower  neutral 
branch  of  the  second  mode  whereas  the 
downstream  branch  point  is  close  to  the  upper 
neutral  branch.  In  the  branch-point  vicinity  (say 
for  the  upstream  point  x14 ),  the  eigenvalues 
behave  as17’22 


a i,2  =  ±  -  xlb  + . . . ,  (24) 

where  ab  =  «,  ( x14  ,a>)  =  a2  (xl4 ,  co) . 

The  receptivity  factor  q  given  by  (16c)  tends  to 
infinity  as  the  resonance  point  tends  to  the  branch 
point,  since  the  scalar  product 

/  dH  \ 

— (x10,«0)A0,B0Uo,  x10^xM.  (25) 

Moreover,  the  eigenvalue  derivative  in  (19)  for 
distributed  receptivity  behaves  as 

7  da„ 

b  =  ~ —  (x10)->°o,  Xi0->%.  (26) 

ax. 

This  singularity  was  analyzed  in  Ref.  [17]  for  the 
case  of  receptivity  to  wall-induced  perturbations. 
Since  the  analysis1  is  easily  extended  to  the  problem 
considered  herein,  we  briefly  outline  key  points  of 
this  extension  for  the  case  of  two-dimensional 
disturbances  ( /3  =  0 ). 

Near  the  branch  point,  one  discrete  mode  (say 
Mode  1  of  the  eigenvalue  a,  (x, ,  co) )  strongly 
interacts  with  another  (say  Mode  2  of  the  eigenvalue 
a2(xt,co) ).  This  interaction  can  be  described  by  the 
two-mode  approximation 

F  (x, ,  y)  =  c,  (x, )  A,  (x, ,  y)  exp  ( A’/lSj )  + 

ci (x, ) A 2 (x, ,  y)  exp( s^'iS2 ) .  (27) 

Far  from  the  resonance  regions,  the  amplitude 
coefficients  c,  2  (x, )  are  solutions  of  the  problem 

dc 

-^  =  ClWu+c2Wn  exp [e~li(S,  -  5, )] ,  (28a) 

dxx 

dc 

=  c2W22  +  cJVn  exp[A'/(.Sj  -  5, )] ,  (28b) 
dxt 

f  (%)  —  cb  ,  c2 (xj  —  c2s .  (28c) 

Flere  cls  and  c2s  are  specified  at  a  certain  initial 
point  x,s ;  the  matrix  elements  are 


B 

’  ’  da  dx , 


B 

1  ’  da  dx, 


where  j  =  1,2  k  =  1, 2  ;  Ht  depends  on  the  mean- 
flow  derivatives  dU  /  5x, ,  dT  /  dx,  and  the  vertical 
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velocity  V0 .  Near  the  branch  point,  the  matrix 
elements  (29)  behave  as 


Wa 


(-1) 


j+k- 1 


4(-ti  -xlb) 


X  — ^  Xti 


(30) 


Then,  the  local  amplitude  coefficient  is 


C„  =  she  1/2  \[xX 


B0,G(x10)F„ 


Bnm,>  (xlb,ab)A( 


fix, 


(35) 


This  leads  to  the  following  singularity  of  the 
amplitude  coefficients 

C  , 

ci,2(*i)  =  7 - ^3— 7j  +  ... ,  x,  — >  xlb ,  (31) 

(*i  ~xib) 

where  constants  C,  2  depend  on  normalization  of 
the  eigenfunctions  A,  2  (x, ,  y) . 

In  the  inner  region  x,  —  xlb  =  0{s212 ) ,  the 
mode  decomposition  (27)  is  replaced  by  the  local 
expansion  with  the  variable  £  =  s~2l2{xl  -xib) 


F  =  ^1/6[c0  (OA0  (y)  +  s~mcm  (^)A01  ( y)  +  ...]x 
exp {ie~'Sb+ie~inab^) ,  (32) 


where  c0  (<Q  is  solution  of  Airy  equation 


d2c0 


A2Cc0=0. 


(33) 


Matching  the  inner  and  outer  solutions  we 
can  establish  an  exchange  rule  in  the  branch¬ 
point  vicinity.21'22  The  local  expansion  (32) 
indicates  that  the  disturbance  amplitude  has  a 
peak  ~£_1/6  in  the  inner  region  =  0(1).  To 
separate  this  effect  from  the  receptivity 
mechanism  we  need  to  establish  a  relationship 
between  the  coefficients  C,  2  (instead  of  the 

amplitude  functions  c12  (x, ) )  and  characteristics 

of  the  incident  acoustic  wave  plus  the  wall 
disturbance. 


Wavy  wall 

In  the  case  of  acoustic  waves  interacting  with 
the  wall  waviness,  the  normal-wave  amplitude  is 
expressed  as 


F„  =  s  he 


(B0,G(x10)Ffl 


Uo  x\b) 


B 


oIIn 

fix, 


(Xi0,  ocQ  )AC 


A„(x1,y)  exp 


*i 

is~{  |  an(xl)dxl  , 
x10 


X,  >  X 


10  5 


(34) 


The  relation  (35)  indicates  that  distributed 
receptivity  has  no  singularity  at  the  spectrum  branch 
point.  In  accordance  with  (31)  the  disturbance  peak  is 
totally  due  to  singular  behavior  of  the  normal  wave 
itself.  If  we  measure  the  disturbance  amplitude  Qn  at 
a  fixed  point  downstream  from  the  branch  point  for 
various  locations  of  the  resonant  point  x10 ,  then  we 
do  not  observe  significant  changes  in  the  distribution 
Q„(x\o )  near  point  x10  =  x16 .  This  conclusion  is 
similar  to  that  formulated  in  Ref.  [17]  for  receptivity 
to  wall-induced  disturbances. 


Local  hump 

If  a  local  hump  is  far  from  the  branch  point,  then 
the  normal-wave  amplitude  is  determined  from  (23). 
As  the  hump  approaches  the  branch  point,  x10  — >  xlb , 
the  receptivity  factor  behaves  as 


C1 


(B0,G(x10)Fa) _ hp(a0) 

B  ^>(x  a  )A  ^ 
ox. 


(36) 


Then,  the  amplitude  coefficient  from  (31)  is 
expressed  as 


C„  = 


£ahAp(a0) 


(B0,G(x10)Fa 


2(x10  Xjj) 


a  (Xio,«o)Ac 
fix, 


•  (37) 


Local  receptivity  increases  proportionally  to 
(x10 -x1A)_1  4  as  the  hump  approaches  the  branch 
point.  This  effect  is  due  to  increasing  of  the  active 
wavenumber  range  of  the  hump  spectrum.  If  the 
hump  center  x10  is  far  from  the  branch  point,  then 

this  range  is  estimated  as  au,+aa-a0  =0(sU2).  If 
Xio -Xu  =  0(e2 /3) ,  then  aw  +  aa -a0  =  0(eU3)  and 
the  amplitude  coefficient  is  Cn  =  0(sahO'6) . 

Matching  of  the  inner  solution  (32)  with  the  outer 
solution  gives  the  normal-wave  amplitude 


F„  =  she 


_,/6  Xp(ah)T{  1/3)  (X 


l4n{xx -xlbf' 


—  I  x 
3 


(B„G(x1A)Fa 

OX , 


■Afl  (x^y)  exp 


is  1  j  an  (Xj  )dx] 


(38) 
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where  s2  3  «  x,  -  xlb  «  1 ;  the  subscript  “ b ” 
denotes  quantities  at  X]  =  xIb ;  the  Gamma- 
function  F(l/3)  =  2.678938...  From  (38)  and 
(3 1 )  we  obtain  the  amplitude  coefficient 


C 


iY/6V(«A)F(l/3)  (B,GFa)4 


/b,^°a\ 

\  ^  A 

(39) 

Summarizing  we  conclude  that  local 
receptivity  strongly  depends  on  the  distance 
between  the  hump  and  the  branch  point.  As 
Xi0  — >  xlb ,  the  amplitude  coefficient  C„ 

—  |x10  —  X\b  |  and  attains  its  maximum 
~£ahe~'16  at  xl0  =  xlb .  If  we  measure  the 
disturbance  amplitude  Q„  at  a  certain  fixed  point 
downstream  from  the  branch  point  for  various 
hump  locations,  then  we  observe  a  local  peak  of 
Q„(xu,)  in  the  vicinity  of  x10  =  xlh .  Again  this 
result  is  consistent  with  that  obtained  in  Ref.  [17] 
for  receptivity  to  wall-induced  disturbances. 

3.  Numerical  results 

Receptivity  calculations  were  conducted  for 
the  boundary  layer  on  a  flat  plate  at  M  =  6 , 
stagnation  temperature  T„  =  390°F  and  the  wall 
temperature  7]*  =  80  F  that  corresponds  to  the 
wall  temperature  ratio  T*/T^=  0.6356.  This 
case  is  a  good  candidate  for  wind  tunnel 
experiments  at  Mach  6.  The  fluid  is  a  perfect  gas 
with  the  specific  heats  ratio  y  =  1.4  and  Prandtl 
number  Pr  =  0.72.  The  viscosity-temperature 
dependency  is  given  by  Sutherland’s  law 

/i{T)  =  (1  +  5)  TV1 ,  (40) 

( T  +  S ) 

where  S  =  110/7],*  for  air  temperature  measured 
in  Kelvin  degrees;  ju  =  ju*  /  fi'e  is  nondimensional 
viscosity.  Hereafter  the  global  length  scale  is 
L"  =  x* ;  the  boundary-layer  scale  is 

8 *  =  A/vf,x*  /  U‘  ,  and  the  Reynolds  number  is 

R=  \lU‘eX  lv*e  ■ 

Effects  of  hypersonic  viscous-inviscid 
interaction  are  neglected;  i.e.,  mean  flow 
characteristics  at  the  upper  boundary-layer  edge 


coincide  with  freestream  characteristics.  The 
boundary-layer  flow  is  approximated  by  the  self¬ 
similar  profiles  U(y),  T(y),  with 

y  =  y*  /  yjv*ex*  l U*  .  The  boundary-layer  thickness  is 

A*  =  16.73^vf*x*  /U*e  and  the  displacement  thickness 

is  A*  =13.64 y]v*x*/U'e  . 

The  second-mode  eigenvalues,  a2(R) ,  are  shown 

in  Fig.  2  for  various  frequency  parameters 

*  *  *2 

F  =  co  veIUe  .  Figures  3  and  4  illustrate  the  phase- 
speed  distributions  c2(R)  =  a>/a2  and  the  branch¬ 
point  loci,  respectively.  Local  maximums 
(minimums)  of  c2(R)  correlate  with  the  upstream 
(downstream)  branch  points.  The  latter  are  close  to 
the  neutral  branches  of  the  second  mode.  The 
upstream  branch  points  are  below  the  real  axis  of 
complex  R  ,  whereas  the  downstream  branch  points 
are  above  this  axis.  In  what  follows,  we  assume  that 
the  branch-point  influence  domains  are  far  enough 
from  the  real  axis.  Then  we  can  use  the  theoretical 
model  of  Section  2  relevant  to  the  simple  discrete 
spectrum. 

Distributions  of  the  pressure  amplitude, 
p2(y)  =  yM2  |A,4(>0| ,  and  the  streamwise  mass-flow 
amplitude, 


My)  = 


A2i  + 


yM2A2i 


(41) 


are  shown  in  Fig.  5  for  77  =  8xl0-5  and  R  =  1000 
(near  the  lower  neutral  branch).  The  eigenfunction  is 
normalized  as  (0)  =  1 .  The  mass-flow  amplitude 
has  a  typical  sharp  maximum  in  the  critical  layer 
located  slightly  below  the  upper  boundary-layer  edge. 


3.1  Waviness-induced  perturbations 

The  pressure  distributions  Re[pui(j/)]  = 

yM2  RefT*^ (j)]  are  shown  in  Fig.  6  for  the 
Reynolds  number  R  =  1000  and  the  wave  numbers 
aw  =  0.001 , 0.05,  0.1.  The  major  portion  of  boundary 
layer  is  weakly  perturbed.  The  first  maximum  of 
pressure  fluctuations  is  located  near  the  upper 
boundary-layer  edge.  The  disturbance  amplitude  in 
the  external  flow  is  essentially  higher  than  that  within 
the  boundary  layer. 

Distributions  of  the  mass  flow  Re[/]y(y)]  are 
shown  in  Fig.  7.  Similar  to  the  case  of  pressure 
disturbances,  the  mass-flow  amplitude  outside  the 
boundary  layer  is  essentially  higher  than  that  within 
the  boundary  layer. 
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3.2  Acoustic  waves 

In  a  supersonic  uniform  flow,  it  is  possible  to 
generate  slow  and  fast  acoustic  waves  with  the 
phase  speeds  ca<  1-1/M  and  ca>  1  +  1/M, 
respectively.  It  is  convenient  to  express  the 
acoustic  wavenumber  as  a  function  of  the  angle 
of  incidence  ©  .  For  slow  waves: 

coM  cos©  co  ,  1  .... 

aa  = - ,  ca  =  —  =  1 - ,  (42) 

Mcos©-l  aa  M  cos© 

0  <  ©  <  ©cr  =  arccos(l/M) , 

For  fast  waves: 

coM  cos©  ,  1  .... 

a= - ,  c  =1  + - ,  (43) 

Mcos©  +  l  M  cos© 

0 < ©  <k 1 2  . 

In  all  figures  presented  hereafter,  ©  is 
measured  in  degrees. 

A  good  insight  into  interactions  of  acoustic 
waves  with  a  hypersonic  boundary  layer  is 
provided  by  the  short-wave  approximation. |Q'20 
Schematic  picture  of  the  disturbance  field  is 
shown  in  Fig.  8  for  slow  waves  with  the  phase 
speed  a(0)  <ca  <Ue-ae,  where  a(y)  is  local 
speed  of  sound.  In  the  lower  wave-guide, 
acoustic  rays  are  reflected  by  the  wall  and 
refracted  near  the  lower  sonic  line, 
U (y)  =  ca  -  a(y) .  In  the  upper  half-space, 
incident  acoustic  waves  are  refracted  near  the 
upper  sonic  line,  U {y)  =  ca  +  a(y) .  In  the  quiet 
zone  between  the  sonic  lines,  the  disturbance 
amplitude  is  relatively  small. 

Slow  waves 

Slow  acoustic  waves  were  calculated  at 
frequency  ^  =  8  x  1 0  and  Reynolds  number 
R  =  1000  .  The  critical  angle  of  incidence  is 
©cr  =  80.4  at  the  Mach  number  M  =  6  .  Figure 
9  shows  contours  of  the  pressure  amplitude 
\p a  (T,  ©)|  •  Figure  10  represents  the  function 
| (jx)|  at  various  angles  of  incidence.  Pressure 

waves  weakly  penetrate  into  the  boundary  layer 
because  its  major  portion  is  covered  by  the  quiet 
zone  (see  Fig.  10).  Due  to  interference  between 
reflected  and  incident  waves,  the  pressure 
amplitude  oscillates  from  0  to  approximately  2 
outside  the  boundary  layer.  This  indicates  that 
the  reflection  coefficient  is  « 1 ;  i.e.,  the 
boundary  layer  is  passive  with  respect  to 
pressure  waves.  A  flat  maximum  of  the  pressure 


amplitude  is  observed  in  the  boundary  layer  at 
©  «  40° .  As  ©  — >  ©„  ,  the  disturbance  wavelength 
decreases,  the  phase  speed  ca  — >  0  .  The  upper  sonic 
line  moves  toward  the  wall  allowing  pressure  waves 
to  penetrate  deeper  into  the  boundary  layer. 

The  contours  of  mass-flow  /j(y,0)  and  the 

functions  |/„(>’)|  are  shown  in  Figs.  11  and  12  for 

various  angles  of  incidence.  As  contrasted  to  the 
pressure  wave,  the  mass-flow  disturbance  has  a  sharp 
peak  near  the  upper  boundary-layer  edge.  Its 
maximum  is  observed  at  ©  «  4(F  .  This  peak  overlaps 
with  the  second-mode  peak  located  in  the  critical- 
layer  (Fig.  5).  Such  an  overlapping  may  increase  the 
scalar  product  (G(jc10  )Fo,B0)  and  enhance 
receptivity. 

Fast  waves 

Characteristics  of  fast  acoustic  waves  were 
calculated  at  the  same  frequency  F  =  8  x  1 0  and 
Reynolds  number  R  =  1000  .  The  pressure  amplitude 
distributions  are  shown  in  Figs.  13,  14.  These  waves 
have  no  sonic  lines  because  their  phase  speed  is 
ca  >1  +  1 1 M  .  Acoustic  rays  easily  penetrate  into  the 
boundary  layer  that  changes  the  pressure  amplitude 
pattern  (compare  Figs.  9  and  13).  The  maximum  of 
\pa(y)\  is  observed  near  the  wall  at  the  angle  of 

incidence  ©  «  50° . 

The  contours  of  mass-flow  |/aCf,©)|  and  the 

functions  |/a(>0|  are  shown  in  Figs.  15  and  16.  The 

mass-flow  amplitude  has  a  peak  near  the  upper 
boundary-layer  edge.  It  is  higher  and  sharper  than  in 
the  case  of  slow  acoustic  waves  (compare  Figs.  12 
and  16). 

Foregoing  numerical  examples  indicate  that 
hypersonic  boundary  layer  is  more  sensitive  to  fast 
acoustic  waves.  This  is  consistent  with  the  acoustic 
field  topology  predicted  by  the  short-wave  theory 
(Fig.  8).  On  the  other  hand,  the  slow- wave  phase 
speed  is  closer  to  the  second-mode  phase  speed  that  is 
more  preferable  for  receptivity.  The  latter  argument  is 
important  for  the  case  of  smooth  surfaces.  However, 
it  is  less  critical  for  receptivity  to  sound  plus  wall 
roughness,  because  the  roughness  shape  can  be 
always  tuned  to  meet  the  synchronization  condition 
(12). 

Small  angles  of  incidence 

Figures  9,  11  (for  slow  waves)  and  Figures  13,  15 
(for  fast  waves)  show  that  the  disturbance  amplitude 
decreases  rapidly  as  the  angle  of  incidence  tends  to 
zero  ( ca  — >1±1  /  M  ).  This  is  due  to  the  fact  that  the 
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reflected  acoustic  wave  is  almost  in  counter- 
phase  with  the  incident  wave  and  they  cancel 
each  other  within  the  boundary  layer.  The  limit 
©  — » 0  needs  special  treatment  because  the 
length  scale  /*  of  the  amplitude  variation  in  the 

vertical  direction  tends  to  infinity.  According  to 
Ref.  [21]  nonparallel  effects  can  change  the 
acoustic  field  substantially  if 

/*/ =  0(AdF'  2)  (in  the  case  considered,  the 
nondimensional  displacement  thickness  is 
Ad  «  14  ).  Using  (42)  we  conclude  that  the  local- 
parallel  approximation  is  not  valid  for 
0  <  0{AdFV2) .  For  disturbances  of  the 

frequency  F  =  8x10^ ,  this  occurs  at  ©  <  7° . 

The  asymptotic  analysis21  showed  that  the 
boundary  layer  growth  (nonparallel  effect) 
causes  diffraction  of  incident  acoustic  waves. 
The  latter  leads  to  decreasing  of  the  acoustic- 
wave  amplitude  at  the  upper  boundary-layer  edge 
in  downstream  direction.  For  acoustic  sources 
fixed  in  space,  the  major  portion  of  the  plate 
surface  is  radiated  by  waves  of  small  ©  .  Then, 
diffraction  and  counter-phase  interference  can 
essentially  reduce  the  sound  amplitude  in  the 
boundary  layer.  This,  in  turn,  reduces  the 
normal-wave  excitation  due  to  sound-roughness 
interaction  far  downstream  from  the  leading 
edge.  Flowever,  in  many  practical  cases,  acoustic 
waves  are  generated  by  convective  sources. 
Their  angles  of  incidence  are  not  small.  For 
example,  Laufer23  showed  that  the  turbulent 
boundary  layer  on  wind-tunnel  walls 
predominantly  radiated  acoustic  waves  with 
©  >  45° . 

These  arguments  indicate  that  receptivity 
experiments  should  cover  a  wide  range  of  angles 
of  incidence.  Special  techniques  need  to  be 
developed  to  simulate  acoustic  field  from 
convective  sources. 


3.3  Distributed  receptivity  on  wavy  surface 
Using  (19)  we  express  the  initial  amplitude  of 
the  second-mode  wave  generated  by  the  wavy 
wall  in  the  resonance  region  as 

\Pl  (*^10 , 0)|  =  £ahGp,dis{X^0}’au)’  (44) 


where  |/>2(x10,0)|  is  pressure  amplitude  on  the 
wall  surface.  The  receptivity  function  is 


Gp,dis  ~  £ 


-1/2,qJ^rM2A24(xw,  0) 


(45) 


The  vector-function  Fa  is  normalized  by  the 
condition:  (pressure  amplitude  of  the  incident  wave  at 
the  upper  boundary-layer  edge)=l;  i.e.  sa  =  p*nc/P* 
at  y  =  ye .  The  nondimensional  amplitude  of  the  wavy 

wall  is  h  =  h *  /  y]v*ex*0  /£/*  ,  where  x*0  is  the  resonance 
point.  The  receptivity  function  Gp  dl!  depends  on  the 
frequency  parameter  F ,  the  Reynolds  number 
R  =  ^julxl  /  v*  and  the  acoustic-wave  incident  angle 
©.  Because  |lm(«2)|  «  |Re(«2)|  for  all  parameters 

under  consideration,  the  resonance  condition  (12)  is 
approximated  as 

«>,  =  «2  (*,„ ,  (o )  -  aa  «  Re[a2  (xl0 ,  co)  -  a  J .  (46) 

Calculations  of  G  dis  were  performed  for  slow 
and  fast  acoustic  waves  of  the  frequency  F  =  8xl0-5 
in  the  Reynolds  number  ranges  relevant  to  the 
second-mode  instability. 

Slow  acoustic  waves 

For  aa  given  by  (42),  the  receptivity  function 
Gpdis(&,R)  is  shown  in  Fig.  17.  This  function  is 

highly  non-uniform  with  respect  to  the  Reynolds 
number.  It  has  two  sharp  maximums  near  the  lower 
and  upper  neutral  branches.  Similar  to  the  case  of 
wall-induced  disturbances,17  these  maximums 
correlate  with  the  spectrum  branch  points  shown  in 
Fig.  4.  The  first  maximum  is  of  the  most  interest, 
since  it  corresponds  to  the  initial  stage  of  instability. 
It  is  observed  at  relatively  high  angles  of  incidence, 
©  «  65' .  Receptivity  to  acoustic  waves  of  small  ©  is 
on  the  order  of  magnitude  lower.  This  is  consistent 
with  the  fact  that  such  acoustic  waves  weakly 
penetrate  into  the  boundary  layer  (see  Section  3.2). 

Fast  acoustic  waves 

For  aa  given  by  (43),  the  receptivity  function 
Gp  dJQ,  R)  is  shown  in  Fig.  18.  Comparing  Figs.  17 

and  18  we  conclude  that  receptivity  to  the  fast 
acoustic  waves  is  higher  than  that  to  the  slow  waves. 
This  agrees  with  the  fact  that  the  fast-wave  amplitude 
is  relatively  large  in  the  boundary  layer  (see  Section 
3.2).  The  first  and  second  maximums  of  Gp  djs  merge 

into  a  dome  covering  the  unstable  range  of  Reynolds 
numbers.  The  receptivity  maximum  is  located  at 
©  ~  45  .  As  the  angle  of  incidence  tends  to  zero,  the 
receptivity  function  decreases  rapidly.  Similar 
features  were  observed  at  other  frequencies  of  the 
range  F  =  (4-10)xl0-5. 
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3.4  Local  receptivity  on  hump 

Using  (23)  we  can  express  the  second-mode 
amplitude  generated  by  acoustic  waves  on  a  local 
hump  as 

\p2{xX0,^  =  sJ^GpM{xl0,co,aa)\p{aw0)\ ,  (47) 

G^^qrM2^  4(*10,0)|,  (48) 

+00 

p{aw)  =  j  g(x')exp(-iawx')dx' .  (49) 


Here  q  is  given  by  (16c);  the  resonance 
wavenumber  aw0  is  approximated  by  (46). 
Similar  to  the  case  of  receptivity  on  a  wavy  wall, 
Gp  1ol  is  a  function  of  the  frequency  parameter 

F ,  the  Reynolds  number  R  =  y]u*exl  /v*  based 
on  the  hump  center  x*0 ,  and  the  acoustic-wave 
incident  angle  © .  The  distribution  of 
(7  ..(©, R)  is  shown  in  Fig.  19  for  slow 

acoustic  waves  and  in  Fig.  20  for  fast  waves. 
Comparing  these  functions  with  the  case  of 
distributed  receptivity  (Figs.  17  and  18)  we 
conclude  that  the  local  receptivity  is  much 
weaker  than  the  distributed  one, 
Gpjoc  ~W  2GPidis.  Nevertheless,  qualitative 
behaviors  of  Gp  loc  and  Gp  dis  are  similar.  This  is 
due  to  the  fact  that  the  ratio, 


does  not  depend  on  characteristics  of  acoustic 
wave  and  roughness. 

In  the  cases  considered,  the  eigenvalue 
derivative  da2  /  dxl  is  a  regular  function  of  xx 
because  the  spectrum  branch  points  are 
essentially  complex  (see  Fig.  4).  The  ratio  (50)  is 
a  smooth  function  of  the  streamwise  coordinate 
x10  or,  what  is  the  same,  the  Reynolds  number 
R  .  However,  for  hypersonic  boundary  layers  on 
a  cooled  wall,17'22  the  branch  points  may  cross 
the  real  axis.  Since  the  eigenvalue  derivative 
tends  to  ±oo  in  the  branch-point  vicinity  (see 
(24)),  the  local  behaviors  of  Gp  djs  and  Gp  loc  are 

different.  In  this  case,  one  needs  to  use  the 
receptivity  model  of  Section  2.3. 


4.  Summary  discussion 


In  this  paper,  the  receptivity  model17  for  wall- 
induced  perturbations  is  extended  to  the  case  of 


acoustic  waves  interacting  with  surface  roughness. 
We  considered  the  second-mode  excitation  on  a  wavy 
wall  (distributed  roughness)  or  a  hump  (local 
roughness  element).  This  problem  includes  the 
following  components: 

(a)  Stationary  disturbances  generated  by  surface 
roughness 

(b)  Acoustic  disturbances 

(c)  Disturbances  resulted  from  the  bilinear  interaction 
of  the  roughness-induced  disturbance  with  the 
acoustic  wave 

It  was  shown  that  the  receptivity  mechanism 
depends  on  the  discrete  spectrum  topology.  If  this 
spectrum  is  simple  (its  branch  points  are  far  from  the 
real  parameter  space),  then  receptivity  reveals  the 
following  features: 

•  For  distributed  roughness  (wavy  wall),  a  strong 
excitation  of  the  normal  waves  occurs  in  the  local 
region  of  the  length  /*  =  R^1:2F  =R1/2<?*,  where 
the  component  (c)  is  in  resonance  with  the 
second-mode  wave.  In  this  region,  the  normal- 
wave  amplitude  is  ~(L*/S*)U2.  The  relative 

length  of  the  resonance  region,  /*/<?*,  increases 
with  the  Reynolds  number.  As  R  — >  oo ,  the 
normal-wave  amplitude  tends  to  infinity.  Its  upper 
bound  is  associated  with  nonlinear  effects  of 
higher-order  approximations. 

•  For  local  roughness,  the  normal-wave  amplitude 
is  finite  even  in  the  parallel  flow  limit.  In  this 
case,  only  a  small  portion  of  roughness  spectrum, 
(a„, -aw0)  =  0(R  ]  2) ,  is  involved  into  the 
resonance.  As  R  — >  oo ,  the  resonant  excitation 
increases  ~  R1  2 ,  whereas  the  active  spectral  band 
is  narrowed  down  ~R~' "2 .  Because  both  trends 
compensate  each  other,  the  normal-wave 
amplitude  tends  to  a  finite  value. 

If  the  branch-point  influence  domain  overlaps 
with  the  real  parameter  space,  then  propagation 
and/or  excitation  of  normal  waves  are  singular. 
Asymptotic  analysis  of  this  singularity  shows: 

•  The  second-mode  excitation  on  a  wavy  wall 
(distributed  roughness)  is  not  singular  at  the 
branch  point.  The  amplitude  peak  is  totally  due  to 
a  singular  behavior  of  the  normal  wave  itself.  If 
we  measure  the  disturbance  amplitude  Q  at  a 

fixed  point  downstream  from  the  branch  point  xxb 
for  various  loci  of  the  resonance  point  x10 ,  then 
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we  do  not  observe  substantial  changes  in  the 
distribution  Q(xw )  near  the  point  xu)  =  xlb . 

•  The  second-mode  excitation  over  a  local 
hump  strongly  depends  on  the  distance 
between  the  hump  center  xl0  and  the  branch 
point  xxb .  As  xw  — >  x]b ,  the  normal-wave 

I  1-1/4 

amplitude  increases  ~  |x10  -  xlb  |  and  attains 
its  maximum  value  at  x10  =  xlb .  If 

we  measure  the  disturbance  amplitude  Q  at  a 
fixed  point  downstream  from  the  branch 
point  for  various  loci  of  the  hump  center,  then 
we  observe  a  local  peak  of  Q(xI0)  in  the 
vicinity  of  x10  =  xlb . 

Numerical  examples  were  presented  for  two- 
dimensional  disturbances  interacting  with  the 
boundary  layer  on  a  flat  plate  with  adiabatic  wall 
at  the  Mach  number  M  =  6 .  In  this  case,  the 
spectrum  branch  points  are  essentially  complex 
that  allows  us  to  calculate  the  receptivity 
coefficients  assuming  that  the  disturbance 
spectrum  is  simple.  Results  of  these  calculations 
are  summarized  as  follows: 

1 .  In  the  boundary  layer,  the  pressure  amplitude 
of  fast  acoustic  waves  is  larger  than  that  of 
slow  waves.  The  amplitude  maximum  is 
observed  at  the  angle  of  incidence  ©  «  40 
for  slow  waves  and  ©  «  45°  for  fast  waves. 
Due  to  counter-phase  interference  between 
the  incident  and  reflected  waves,  the 
disturbance  amplitude  decreases  rapidly  as 
©^0. 

2.  On  wavy-wall,  the  receptivity  coefficient  is 
highly  non-uniform  with  respect  to  the 

Reynolds  number  R  =  yju*exl  /  v*  . 

•  For  slow  acoustic  waves,  the  receptivity 
coefficient  has  two  maximums,  which  are 
close  to  the  lower  and  upper  neutral  branches. 
Loci  of  these  maximums  correlate  with  the 
spectrum  branch  points.  The  boundary  layer 
is  most  receptive  to  acoustic  waves  of 
relatively  high  angles  of  incidence:  ©  «  65° 
near  the  first  maximum  and  ©  «  40°  near  the 
second  maximum. 

•  For  fast  acoustic  waves,  the  receptivity 
maximums  merge  into  a  dome  covering  the 
instability  region.  The  boundary  layer  is  most 
receptive  to  acoustic  waves  of  ©  «  45° . 


•  The  receptivity  coefficients  for  both  slow  and  fast 
acoustic  waves  of  small  angles  of  incidence  are  on 
the  order  of  magnitude  smaller  of  their  maximal 
values. 

3.  On  local  hump,  the  receptivity  coefficient  has  the 
same  qualitative  behavior  as  in  the  wavy-wall 
case.  Flowever  its  quantity  is  essentially  smaller: 

GpJoc  *10 -2GPtdis. 
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Fig.  1  Flow  scheme. 
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Fig.  2  Second-mode  eigenvalues  a2  as  a  function  of 
Reynolds  number  at  various  frequencies. 
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Fig.  3  Second-mode  phase  speed  as  a  function  of 
Reynolds  number  at  various  frequencies. 


Fig.  6  Distributions  Re[plv(y)]  for  wall-induced 
perturbation  at  various  wavenumbers  a  =  aw ; 
=  1000. 


Fig.  4  Blanch  points  in  complex  R  -plane  at  pjg  7  Mass-flux  distributions  Re(/W(>’))  for  wall- 

various  frequencies.  .  ,  ,  ,  . 

induced  perturbation  at  various  wavenumbers  a  =  aw  ; 


R  =  1000  . 


Fig.  5  Second-mode  eifenfunctions  for  pressure 

p2(y)  and  streamwise  mass  flux  f2{y)  near  the  Fig'  8  Scheme  of  acoustic  fleld  in  the  boundary  layer, 
lower  neutral  curve:  R  =  1000  ,  F  =  8.x  10  3 . 
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Fig.  9  Contours  |/>o(y,©)|  =  const  for  slow  waves; 
£  =  1000,  F  =  8.xl0^5 . 


Fig.  12  Distributions  |/„(j,@)|  at  various  incident 
angles  for  slow  waves  R  =  1000  ,  F  =  8.x  10"5 . 


Fig.  10  Distributions  |/?a(j,©)|  at  various  incident 
angles  for  slow  waves,  R  =  1000  ,  F  =  8.x  10"5 . 


Fig.  13  Contours  \pa(y,®)\  =  const  for  fast  waves; 
£  =  1000,  £  =  8.xl0~5. 


Fig.  11  Contours  of  \fa(y,&)\  =  const  for  slow 
waves;  £  =  1000,  £  =  8.xl0~5. 


Fig.  14  Distributions  |/>a(y,©)|  at  various  incident 
angles  for  fast  waves,  £  =  1000  ,  F  =  8.x  10  5 . 
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R  =  1000  ,  F  =  8.x  10-5 .  Fig.  18  Distributed  receptivity  function  for  fast  acoustic 

waves;  F  =  8.xl0~5 . 
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